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I,  Introduction 

y  This  final  report  covers  the  most  recent  period  of  the  Office  of  Naval 
Research  Contract  NgT  N00014-75-C-0245 .^titled' ^DYNAMICAL  PROPERTIES  OF 
DISORDERED  SYSTEMS  .  During  this  period,  research  has  been  completed  on 
a  variety  of  research  problems  contained  in/iONR  Research  Proposal.  The 
principal  areas  of  research  can  be  summarized  under  the  following  categories: 

1)  The  fractal  interpretation  of  amorphous  structures 

2)  Non-linear  conductivity  of  low  dimensional  materials, 

3)  The  effect  of  anodization  on  the  superconducting  transition  temperature, 

4)  The  time  dependence  of  the  remanent  magnetization  of  spin  glasses 

5)  The  experimental  opportunities  for  physics  at  high  magnetic  fields  *  f'. _ 

In  the  next  Section,  we  outline  the  research  accomplishments  in  each  of 
these  areas.  The  detailed  references  to  published  research  supported  in 
these  areas  by  the  Office  of  Naval  Research  is  contained  in  Section  III. 

II.  Research  accomplishments 

1)  The  fractal  interpretation  of  amorphous  structures 

During  the  period  1982-83,  we  began  an  intensive  investigation  of  the 
dynamical  properties  of  physical  systems  which  were  fractal  in  geometry.  We 
showed  that  the  density  of  states  on  a  fractal  can  be  calculated  by  taking 
into  account  the  scaling  properties  of  both  the  volume  and  the  connectivity. 

We  used  a  Green's  function  method  developed  elsewhere,  utlizing  a  relationship 
to  the  diffusion  problem.  We  found  that  proper  mode  counting=required  a 
reciprocal  space  with  a  new  intrinsic  fracton  dimensionality  d  =  2d/(2  +  0). 

Here,  d  is  the  effective  dimensionality,  and  0  the  exponent  giving  the  dependence 
of  the  diffusion  constant  on  distance.  For  example,  we  find  for  percolation 
clusters  d  =  4/3  within  the  numerical  accuracy  available,  independent  of 
Euclidean  dimensionality  d.  Crossover  to  normal  behavior  at  low  frequencies 
was  discussed  for  finite  fractals  and  for  percolation  above  the  percolation 
threshold  pc>  Relevance  to  experimental  results  on  proteins  was  also  discussed. 

In  a  follow-up  piece  of  research,  we  calculated  the  density  of  states 
for  lattice  vibrations  on  a  fractal  with  careful  attention  paid  to  the 
normalization  condition.  It  was  found  that  at  the  cross-over  between  Debye-type 
excitations  (long  wavelength)  and  "fracton"  excitations  (short  length  scale) 
the  density  of  states  is  discontinuous.  The  size  of  the  discontinuity  was 
related  to  the  ratio  of  the  fracton  dimensionality  to  the  Euclidean  dimensionality. 
An  application  was  made  to  percolating  structures.  A  set  of  missing  modes  was 
identified  which  may  be  the  origin  of  the  two-level  systems  hypothesized  for 
amorphous  structures.  The  specific  heat  of  epoxy-resin  exhibits  a  cross-over 
from  a  Debye-type  region  (T<8K)  to  a  region  (8-50K)  where  the  vibrational  density 
of  states  depends  linearly  on  the  frequency.  Over  the  same  frequency  regime, 
the  thermal  conductivity  exhibits  an  effective  phonon  mean  free  path  of  the 
order  of  (or  less  then)  a  lattice  constant.  We  interpreted  this  behavior  in 
terms  of  quantized  fractons,  with  an  energy  range  8-50K,  and  we  suggested  that  these 
fracton  states  are  localized.  This  was  consistent  with  the  usual  interpretation 
of  a  precipitous  drop  in  the  phonon  mean  free  path  at  the  cross-over  energy 
of  8K.  Analogous  behavior  was  argued  for  the  thermal  properties  of  glasses  which 
exhibit  a  similar  structure  in  the  thermal  conductivity.  Recent  neutron- 
irradiated  quartz  experiments  tend  to  confirm  this  interpretation. 


2)  Non-linear  conductivity  of  low  dimensional  materials 

The  frequency  dependence  of  the  electrical  conductivity  was  examined  for 
a  classical  hopping  model  of  a  random  one-dimensional  system  in  the  presence 
of  a  superposed  static  electric  field.  The  effect  of  the  field  was  taken  as 
a  constant  bias  for  the  left-right  jump  rates.  A  general  expression  was 
derived  for  the  mean  velocity  and  frequency-dependent  conductivity.  Explicit 
evaluation  of  these  equations  was  given  for  correlated  and  uncorrelated 
hopping  rates:  (1)  in  general  for  high  frequencies  and  (2)  to  lowest  order 
in  the  disorder  for  all  frequencies.  In  the  latter  case,  an  initial  decrease 
in  the  frequency  dependent  conductivity  for  very  small  frequencies  was  found, 
with  the  real  part  of  the  conductivity  varying  as  the  square  of  the  frequency 
while  the  imaginary  pary  vary led  as  the  cube  of  the  frequency.  For  larger 
frequencies,  the  conductivity  crossed  over  to  the  form  previously  calculated 
by  Alexander  and  Orbach  (proportional  to  the  square  root  of  the  frequency). 
The  cross-over  frequency  increases  with  the  bias.  In  addition,  the  variance 
in  the  autocorrelation  function  was  calculated  in  the  long-time  limit,  for 
weak  disorder  in  the  symmetric  case.  It  was  shown  that  fluctuations  do  not 
significantly  affect  the  determination  of  this  quantity  under  these 
conditions. 

3)  The  effect  of  anodization  on  the  superconducting  transition 
temperature 

We  measured  the  superconducting  transition  temperature  of  anodized 
aluminum  films  of  grain  sizes  ranging  from  less  than  100  to  3,000  angstroms. 
The  transition  temperature  is  1.8K  for  films  of  grain  size  100  angstroms,  and 
decreased  monotonically  with  increasing  grain  size  to  1.2K  for  3,000  angstrom 
grains.  The  effect  depended  only  on  the  volume  of  the  grains.  We  concluded 
that  the  phonon  softening  model  appeared  to  break  down  for  thin  aluminum 
grains. 

4)  The  time  dependence  of  the  remanent  magnetization  of  spin  glasses 

The  time  decay  of  the  thermoremanent  magnetization  (M_„m)  was  measured 

in  1.0%  Cu:Mn  and  2.6%  Ag:Mn  spin-glasses.  It  was  shown  tnar  was 

neither  an  algebraic  nor  a  logarithmic  function  of  time.  InsteaaY  it  was 
found  that  H-RM  can  be  characterized  by  a  "streched"  exponential: 

*  Mgexp[--C(t/tQ)  n/(l-n)].  The  prefactor  (Mq)  and  the  time-stretch 
exponent  (n)  are  temperature  dependent,  whereas  the  exponential  factor  (C) 
and  the  relaxation  rate  (1/t^)  can  be  chosen  to  be  independent  of  temperature 
throughout  the  spin  glass  region.  In  addition,  we  found  Mq  to  depend 
linearly  on  the  cooling  field  H,  whereas  n,  C,  and  1/tQ  are  independent  of  H 
for  H  less  than  30  Oersteds.  Similar  time  dependences  appear  in  the 
disorder-diffusion  theory  of  Grassberger  and  Procaccia  and  the 
cooperative-relaxation  theory  of  Ngai,  but  neither  theory  in  its  present  form 
was  found  to  be  directly  applicable  to  spin  glasses. 

5)  The  experimental  opportunities  for  physics  at  high  magnetic  fields 

A  brief  list  of  current  areas  of  research  in  high  magnetic  field  physics 
was  investigated.  Three  distinct  areas  were  examined  in  detail:  (1)  density 
of  states  for  vibrational  states  on  fractals  ("f ractons") ,  (2)  the 
thermodynamic  properties  of  exchange  enhanced  quasi-magnetic  systems,  and  (3) 
p-state  pairing  in  thin  film  or  layered  superconductors. 
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R6sum6.  —  Nous  calculons  la  denstte  detats  sur  un  fractal  en  tenant  compte  des  proprietes  dechelles 
pour  le  volume  et  la  connectivity  Nous  utilisons  la  methode  de  Green  developpee  par  ailleurs, 
qui  utilise  une  relation  au  probleme  de  diffusion.  Nous  avons  trouve  que  pour  compter  les  modes 
correctement,  on  doit  avoir  un  espace  reciproque  avec  une  nouvelle  dimensionnalite  de  fracton 

4  =  2  4/(2  +  8) . 

Ici  4  est  la  dimension  effective,  et  J  est  l'exposant  que  caracterise  la  variation  de  la  constante  de 
diffusion  avec  la  distance.  Par  exemple.  nous  trouvons  pour  les  amas  de  percolation,  4  =  43,  quelle 
que  soit  la  dimensionnaliti  Euclidienne  d  i  la  precision  numerique  disponible.  Nous  discutons  le 
«  crossover  »  vers  un  comportement  normal  aux  basses  frequences  pour  des  fractals  finis,  et  pour 
la  percolation  au-dessus  du  seuil  de  percolation  pc.  Nous  examinons  aussi  la  pertinence  de  nos 
predictions  en  les  confrontant  k  des  resultats  experimental  sur  les  proteines. 

Abstract  —  The  density  of  states  on  a  fractal  is  calculated  taking  into  account  the  scaling  properties 
of  both  the  volutqe  and  the  connectivity.  We  use  a  Green's  function  method  developed  elsewhere 
which  utilizes  a  relationship  to  the  diffusion  problem.  It  is  found  that  proper  mode  counting  requires 
a  reciprocal  space  with  new  intrinsic  fracton  dimensionality  J  =  2  4/(2  +  J).  Here,  4  is  the  effective 
dimensionality,  and  J  the  exponent  giving  the  dependence  of  the  diffusion  constant  on  distance. 
For  example,  we  find  for  percolation  clusters  4  »  4/3  within  the  numerical  accuracy  available, 
independent  of  the  Euclidean  dimensionality  4.  Crossover  to  normal  behaviour  at  low  frequencies 
is  discussed  for  finite  fractals  and  for  percolation  above  the  percolation  threshold  pc.  Relevance 
to  experimental  results  on  proteins  is  also  discussed. 


(•)  Permanent  address  :  Department  of  Physics,  University  of  California,  Los  Angeles,  California. 
90024,  U.S.A. 
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The  purpose  of  this  note  is  to  calculate  the  density  of  states  on  fractals.  We  do  this  using  a 
Green's  function  technique  [1],  which  avoids  the  use  of  boundary  conditions  and  wave  vector 
counting.  We  show  that  the  usual  relationship  between  the  density  of  states  and  the  Euclidean 
dimensionality  {d  I  does  not  apply.  The  density  of  states  cannot  be  described  in  terms  of  an  ano¬ 
malous  dimensionality  (d)  alone  but  requires  an  additional  index  describing  the  internal  struc¬ 
ture.  We  determine  a  fracton  dimensionality  (d)  of  the  relevant  reciprocal  space  which  assures 
proper  mode  counting  for  a  fractal  in  terms  of  the  index  governing  diffusion  (J)  and  the  ano¬ 
malous  or  fractal  dimension  ( d ).  The  density  of  states  for  free  particles  and  for  lattice  vibrations 
is  determined.  We  apply  these  results  to  polymer  chains,  the  triangular  Sierpinski  gasket  and 
to  percolation  networks. 

We  were  directly  motivated  by  the  recent  work  of  Stapleton  et  al.  [2]  who  found  an  anoma¬ 
lous  temperature  dependence  for  the  ESR  spin-lattice  relaxation  time  of  iron  in  several  proteins. 
They  interpreted  their  measurements  in  terms  of  an  anomalous  vibrational  density  of  states 
arising  from  a  suggested  fractal  structure  for  the  proteins.  However,  they  included  only  the 
anomalous  dimensionality  d  in  their  analysis. 

The  density  of  states  on  fractals  should  also  be  of  interest  in  other  situations.  Examples  would 
be  the  role  of  geometrical  disorder  in  amorphous  systems,  and  the  specific  heat  of  fractal-like 
systems  in  the  intermediate  temperature  range  (see  below). 

1.  Method  of  calculation.  —  We  consider  problems  defined  on  a  fractal  so  that  both  the 
available  volume  and  the  connectivity  are  determined  by  the  fractal  geometry.  An  explicit  way 
of  realizing  a  fractal  model  is  to  consider  solutions  of  the  relevant  equations  on  a  network  of 
wires  connected  in  a  suitable  geometry  [3].  We  take  the  network  to  be  homogeneous.  Our  results 
will  apply  for  times  or  frequencies  such  that  the  associated  distances  are  much  larger  than  the 
size  of  the  individual  bonds  making  up  the  network,  but  much  smaller  than  the  total  size  of  the 
network.  Formally,  we  are  assuming  that  one  can  define  some  local  geometry  on  the  fractal  in 
which  the  Laplacian  has  its  usual  form  (i.e.  is  equivalent  to  a  local  q 2  expansion ).  This  is  defined 
by  the  wires  for  a  network  model.  While  this  may  not  be  possible  for  all  conceivable  fractals, 
we  believe  it  is  certainly  possible  for  most  cases  in  which  one  would  be  interested.  One  finds 
that  the  formal  relationship  between  problems  initially  described  by  Laplacians  (or  by  equi¬ 
valent  finite  difference  equations  (4j)  is  maintained. 

The  structure  of  the  diffusion  equation  is  such  that  it  can  be  mapped  onto  a  master  equa¬ 
tion.  which  in  turn  has  the  same  form  as  the  free  particle  Schrodinger  equation  and  the  equa¬ 
tion  of  motion  for  mechanical  vibrations  (see  Sec.  3  of  Ref.  [3],  and  Ref.  [5]).  This  will  enable  us 
to  map  the  eigenvalue  density  of  states  for  the  quantum  vibrational  problem  onto  the  eigen¬ 
value  density  of  states  of  the  diffusion  problem.  The  latter  can  be  obtained  from  the  single  site 
Green's  function  for  the  diffusion  problem  [1]  : 

N(e)  =  -  ^  Im  <  P0(- £  +  i0‘)  >  (1) 

where  P0{e)  is  the  Laplace  transform  of  P0(rL  the  autocorrelation  function,  with  e  the  spectral 
parameter.  In  physical  terms,  if  a  particle  is  initially  at  the  origin  at  time  t  =  0,  the  probability 
of  finding  it  there  at  time  t  is  given  by  P0(r ).  We  now  calculate  this  quantity  directly. 

On  a  fractal  one  expects,  in  general,  anomalous  diffusion  [6J.  We  write. 

<  rJ(f)  >  x  t2,2+*\  (2) 

where  d  is  an  index  which  depends  on  the  geometry  (i.e.  it  is  the  diffusion  constant  scaling  expo¬ 
nent).  In  general  one  expects  6  *  0.  The  total  volume  available  on  the  fractal  within  the  dif¬ 
fusion  distance,  is 

F(t)  x  <  r2(r)>?:,  (3) 
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where  d  is  the  anomalous  or  fractal  dimensionality.  Thus. 

<  P0(f)>  x  [TUI]'1  x  (4) 

We  have  thus  assumed  that  the  diffusion  length  is  given  by  a  power  law  (Eq.  (2))  and  that 
it  is  the  only  relevant  length  scale  in  the  problem.  Using  equation  (1).  one  is  immediately  led 
to  the  eigenvalue  density  of  states  for  the  diffusion  problem  : 

N{ e)  x  r* ,  (5) 

where 

x  =  [1(2  +  J)]  -  1  .  (6) 

As  shown  explicitly  in  reference  (1],  the  spectral  parameter  e  can  be  related  to  the  energy  eigen¬ 
values  of  the  vibrational  problem  by  replacing  e  by  to2  and  multiplying  by  to  : 

N(u>)  x  to'  ,  (7) 

where 

p  =  2  x  +  1  =  [2  5/(2  +  £)]  -  1  .  (8) 

We  refer  to  thequantized  vibrational  states  on  a  fractal  as  fractons. 

For  positive  6,  x  and  p  are  always  smaller  than  the  anomalous  dimension  result  (x  =  (d  2)-  1, 
p  -  d  -  1,  as  suggested  in  Ref.  [2]).  The  reason  for  the  difference  is  that  the  scale  dependence 
of  the  elastic  constants  is  also  anomalous  (i.e.  it  depends  on  6).  It  is  not  sufficient  to  consider 
only  the  mass  scaling  described  by  d. 

Our  results,  equations  (5)-(8).  can  be  described  by  a  «  mode  counting  »  reciprocal  space  of 
effective  dimensionality 

d  =  2  d/( 2  +  J) .  (9) 

We  shall  cal!  this  the  fracton  dimensionality.  For  a  standard  Laplacian  (d  -»  q2)  expansion, 
this  is  the  relevant  dimensionality.  It  determines  the  relevant  Hilbert  jpace  for  Laplacian  equa¬ 
tions  on  the  fractal.  We  emphasize  that  the  fracton  dimensionality  (d )  is  an  intrinsic  property 
of  the  fractal  geometry.  It  differs  from  the  mass  scaling  exponent,  or  the  fractal  dimensiona¬ 
lity  (dl  and  from  the  diffusion  constant  scaling  exponent  (<H  in  being  independent  of  the  manner 
in  which  the  fractal  is  embedded  in  an  external  space  (of  Euclidean  dimension  d).  The  d  depen¬ 
dence  of  d  and  S  cancels  out  in  equation  (9). 

We  note  that  this  seems  to  be  the  natural  extension  of  a  gradient  expansion  to  fractals.  The 
obvious  alternative  of  expanding  in  undistorted  real  space  gradients  leads  to  a  singular  expan¬ 
sion  (e  -  q2*7)  and  a  reciprocal  space  of  dimensionality  d  (using  Eq.  (2)). 

2.  Some  examples.  —  It  is  useful  to  consider  some  examples. 

a)  Consider  a  one  dimensional  chain  whose  configuration  is  described  by  a  random  walk 
or  self-avoiding  walk.  One  has 

N’  x  r.  (10) 

where  N  is  the  length  (number  of  units)  of  the  chain.  Thus. 

d  =  v-'  .  (11) 


For  diffusion  along  the  chain. 


N!(i|xi, 
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and  from  equations  (10)  and  (11). 

<  r2(r)  >  x  <  JV2v(r)  >  x  t '  . 


so  that 


8  =  2(d  -  1), 
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(12) 

(13) 


leading  to  (Eqs.  (6),  (8),  and  (9)), 

x  =  1/2,  p  =  0,  and  d  =  1  .  (14) 

The  fracton  dimension  (d )  for  this  problem  is  identically  one,  and  is  independent  of  v. 

b)  Consider  a  triangular  (d  =  2)  Sierpinski  gasket.  One  has  [3.  7] 

8  =  [(In  5)/ (In  2)]  -  2  *  0.322 ,  (15a) 

d  =  (In3)/(ln2)  *  1.585,  (156) 

but  the  fracton  dimension  is  (Eq.  (9)), 

d  =*  2(ln  3)/ (In  5)  1.365  .  (16) 

Thus,  the  gasket  is  somehow  more  one  dimensional  than  suggested  by  d.  For  the  eigenvalue 

density  of  states,  one  has  the  spectral  parameter  exponent  (Eq.  (6)), 

x  =  [(In  3)/ (In  5)]  -  I  %  -  0.317 ,  (17 a) 

and  the  fracton  oj  exponent  (Eq.  (8)), 

p  =  [2(ln  3)/ (In  5)]  -  1  %  1.365  .  (176) 

c)  Consider  a  critical  percolation  network.  Taking  the  infinite  cluster  at  the  critical  percola¬ 
tion  concentration  pc,  one  has  from  straight  forward  scaling  considerations  [8], 

d  =  d  -  (0/v),  (18) 

where  d  is  the  Euclidean  dimensionality  of  the  external  space,  and  (3,  6] 

J=(f-/?)/v,  (19) 

where  f  is  the  conductivity  exponent.  Thus, 

d  =  2(dv  -  P)i(t  -  p  +  2  v) .  (20) 


Table  I.  —  The  fracton  dimensionality  d  defined  in  equation  (9)  and  the  eigenvalue  density  of 
states  indices  x  (Eq.  (6))  and  p  (Eq.  (8))  as  a  function  of  the  Euclidean  dimensionality  of  the  per¬ 
colation  problem  (d).  The  fractal  dimensionality  d  (Eq.  (18))  and  8  (Eq.  (19))  were  computed  from 
the  numerical  values  in  reference [8]. 


d 

d 

X 

P 

d 

8 

2 

1.36 

-  0.32 

0.36 

1.9 

0.80 

3 

1.42 

-  0.29 

0.42 

2.5 

1.55 

4 

1.39 

-  0.30 

0.39 

3.3 

2.71 

5 

1.44 

-  0.28 

0.44 

3.8 

3.3 

X 

4/3 

-  1/3 

1/3 

4 

4 
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We  give  results  for  d  in  table  I  as  a  function  of  the  Euclidean  dimensionality  (d)  for  Stauffer's 
values  [8]  of  the  indices  /,  f},  and  v  forjhe  percolation  problem.  One  notes  the  very  weak  depen¬ 
dence  of  the  fracton  dimensionality  (d)  on  the  Euclidean  dimensionality  (d)  of  the  percolation 
problem.  This  is  in  sharp  contrast  to  the  behaviour  of  d  and  8.  The  table  certainly  suggests  the 
conjecture  that,  for  percolation  on  the  infinite  cluster, 

1=4/3,  (21) 

independent  of  d  (see  note  added  in  proof). 

3.  Crossover  and  finite  size  effects.  —  The  Laplace  transform  of  equation  (2)  yields  a  rela¬ 
tion  between  length  scale  and  the  Laplace  transform  spectral  parameter  t.  We  are  able  to  find 
a  characteristic  diffusion  length  scale  appropriate  to  the  spectral  parameter  £,  which  varies 
as  [3] 

;.toc  £-,/U+J|.  (22a) 

Mapping  onto  the  vibrational  problem  as  before, 

kux.<j)-1:a*h .  (22 b) 

Fractal  behaviour  is  found  for  length  scales  less  than  the  size  (L)  of  a  fractal  object.  For  larger 
length  scales,  the  solutions  are,  in  essence,  uniform  over  the  fractal  and  the  density  of  states  is 
determined  by  the  boundary  conditions.  Crossover  to  fractal  behaviour  occurs  when  /.t  $  L, 
or  for  energies 

«c.o L-,2+J\  (23) 

where  the  subscript  c.o.  means  crossover.  For  e  >  ec  o  (or  to  >  coc  o  ),  the  eigenvalue  density 
of  states  is  then  given  by  equation  (5)  (or  Eq.  (7)).  For  a  percolation  network  above  pc.  one  pre¬ 
dicts  normal  d-dimensional  low  frequency  vibrational  density  of  states  (p  =  d  -  1)  crossing 
over  to  fracton  behaviour  for  vibrational  frequencies  above  (Eq.  (23)) 

ojc.0.(p)  cc  <Jp,2  +  a''2  (24) 


where  <JP  is  the  percolation  correlation  length. 

4.  Discussion.  —  We  conclude  with  some  remarks  concerning  the  relevance  of  our  results 
to  experiment.  We  have  analysed  the  eigenvalue  density  of  states  for  diffusion  on  a  fractal  geo¬ 
metry,  and  mapped  our  results  onto  the  vibrational  eigenvalue  density  of  states  for  systems 
with  the  same  geometrical  structure.  For  the  vibration  problem,  this  assumes  that  both  the 
elastic  constants  and  the  inertial  mass  are  appropriate  to  a  free  fractal.  It  is  important  to  empha¬ 
size  that  this  is  not  necessarily  a  proper  physical  description.  It  is,  in  fact,  hard  to  think  of  a 
situation  where  the  vibration  spectrum  of  a  polymer  would  be  adequately  described  by  our 
free  fracton  model.  For  gels,  the  mass  density  is  always  dominated  by  the  solvent  and  there¬ 
fore  scales  with  the  Euclidean  dimensionality  of  the  external  space  ( d ).  Any  fractal  anomalies 
would  only  reflect  the  scaling  properties  of  the  elastic  constants.  It  is  conceivable  that  something 
similar  occurs  for  proteins  where  one  would  also  expect  an  essentially  uniform  density.  This 
might  be  relevant  to  the  interpretation  of  the  results  of  Stapleton  el  al.  [2]  and  to  other  measu¬ 
rements  in  which  the  vibrational  density  of  states  of  proteins  is  important  [9], 
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Office  of  Naval  Research.  We  wish  to  acknowledge  helpful  correspondence  with  H.  J.  Stapleton, 
and  very  useful  conversations  with  T.  C.  Lubensky.  R.  Rammal,  G.  Toulouse,  and  J.  Vanni- 
menus  (see  note  added  in  proof). 

Note  added  in  proof.  —  Professor  T.  C.  Lubensky  has  noted  (private  communication)  that 
some  intriguing  consequences  follow  if  one  takes  seriously  the  results  of  the  table  for  the  per¬ 
colation  problem.  He  notes  that  if  one  assumes  that  d  =  4 .3,  independent  of  d,  one  can  use 
equation  (20)  to  generate  an  expression  for  the  exponent  t  in  terms  of  v  and  0  : 

r  =  (1/2) [v(3 d-  4)  -  ft]. 

Two  consequences  are  worth  noting. 

1)  Using  the  values  of  /?  and  v  following  from  the  den  Nijs  conjecture  at  d=2  (M.  P.  M.  den 
Nijs,  Physica  (Utrecht)  A  95  (1979)  449)  :  /J  =  5/36  and  v  =  4/3.  one  finds  t  =  91/72  =  1.264. 

This  differs  from  t  =  1.1  presented  by  S.  Kirkpatrick  (La  matiere  mal  condensee.  Ed.  by  R.  Ba- 
lian,  R.  Maynard  and  G.  Toulouse  (North-Holland  publishing  company,  Amsterdam)  1979, 
p.  321),  and  is  not  in  accord  with  the  relation  t  =  v  (A.  K.  Sarychev  and  A.  P.  Vinogradoff  J. 
Phys.  C  14  (1981)  L-487).  A  very  recent  finite  size  scaling  (simulations  on  finite  size  strips)  of 
B.  Derrida  and  J.  Vannimenus  (submitted  for  publication,  1982)  finds  t  =  1.28,  in  close  accord 
with  the  above  consequence  (r  =  1.264)  of  setting  d  =  4,3. 

2)  The  links  and  node  model  led  A.  S.  Skal  and  B.  I.  Shklovskii  [Fiz.  T ekh.  Poluprov.  8  (1974) 
1582  (Sov.  Phys.  :  Semicond.  8  (1975)  1029)]  to  define  t  -  (d  -  2)  v  +  {.  Setting  d  =  4/3,  solving 
for  r,  and  using  appropriate  scaling  relationships,  leads  to  the  expression,  valid  for  all  d, 

C  =  (1/2)  (p  +  y). 

G.  S.  Grest  and  M.  J.  Stephen  (Phys.  Rev.  Lett.  38  (1977)  567)  and  C.  Dasgupta,  A.  B.  Harris 
and  T.  C.  Lubensky  (Phys.  Rev.  B  17  (1978)  1375)  show  that  =  1  +  0(e2)  where  e  -  6  -  d. 
If  we  use  perturbation  theory  results  for  /J  and  y  (R.  G.  Priest  and  T.  C.  Lubensky,  Phvs.  Rev. 
B  13  (1976)  4159 ;  B  14  (1976)  5125 ;  D.  J.  Amit,  J.  Phys.  A  9  (1976)  1441)  : 

>5=1  —  (1/7)6  -  (6U73  32  22)e2  +  - 

y  =  1  +  (1/7)  e  +  (565/73  32  22)  e2  +  - 

we  find 

Cp«„  =  1  +  CM9 

as  required  to  0(e2).  For  4  =  2,  t^  =  £  =  1.33,  obviously  outside  the  range  of  «  safe  »  conver¬ 
gence  (e  =  4),  but  remarkably  close  to  the  conjecture  (d  =  4/3)  value  of  t  =  1.264. 

Finally,  from  another  perspective,  R.  Rammal  and  Angles  d'Aurioc  (submitted  for  publica¬ 
tion,  1982)  have  found  values  for  d  and  5  for  the  Sierpinski  sponge  in  d  dimensions.  They  find 

d  =  [In  (d  +  l)]/ln  2,  5  =  { [In  (d  +  3)]/ln  2  }  -  2 

so  that 

1  =  2[ln  (d  +  l)]/ln  (d  +  3) . 

These  results  agree  with  ours  for  the  case  we  considered,  d  =  2.  In  addition,  the  asymptotic 
d  x>  limit  for  d  is  2,  quite  different  from  the  d  -*  x  limit  for  the  percolation  problem,  4/3. 
The  ratio,  d/d  =  (In  4)/[ln  (d  +  3)]  is  always  less  than  unity  for  d  ^  Z  exhibiting  the  same 
trend  as  we  have  exhibited  for  d  -  2. 
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Fracton  interpretation  of  vibrational  properties  of  cross-linked  polymers, 
glasses,  and  irradiated  quartz 
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The  density  of  states  for  thermal  vibrations  on  a  fractal  is  calculated  with  careful  attention  paid 
to  the  normalization  condition.  It  is  found  that  at  the  crossover  between  Debye-type  excitations 
(long  wavelength)  and  “fracton"  excitations  (short-length  scale)  the  density  of  states  is  discontinu¬ 
ous.  The  size  of  the  discontinuity  is  related  to  the  ratio  of  the  fracton  dimensionality  to  the  Euclide¬ 
an  dimensionality.  Application  is  made  to  percolating  structures.  A  set  of  missing  modes  is  identi¬ 
fied  which  may  be  the  origin  of  the  two-level  systems  hypothesized  for  amorphous  structures.  The 
specific  heat  of  epoxy  resin  exhibits  a  crossover  from  a  Debye-type  region  I  T  <  8  K)  to  a  region 
(8—50  K)  where  the  vibrational  density  of  states  depends  linearly  on  the  frequency.  Over  the  same 
frequency  regime,  the  thermal  conductivity  exhibits  an  effective  phonon  mean  free  path  of  the  order 
of  (or  less  than)  a  lattice  constant.  We  interpret  this  behavior  in  terms  of  quantized  fractons,  with 
an  energy  range  8—50  K,  and  we  suggest  that  these  fracton  states  are  localized.  This  is  consistent 
with  the  usual  interpretation  of  a  precipitous  drop  in  the  phonon  mean  free  path  at  the  crossover  en¬ 
ergy  of  8  K  Analogous  behavior  is  argued  for  the  thermal  properties  of  glasses  which  exhibit  a 
similar  structure  in  the  thermal  conductivity.  Recent  neutron-irradiated  quartz  experiments  tend  to 
confirm  this  interpretation. 


I.  INTRODUCTION 

Two  of  us  have  suggested  that  the  concept  of  fractals1 
can  be  applied  to  the  vibrational  properties  of  macro¬ 
molecules  and  have  derived  the  vibrational  density  of 
states  for  fractal  structures.2  We  found  that  the  usual 
Debye-type  density  of  states  crosses  over  to  a  "fracton" 
density  of  states  for  length  scales  less  than  some  charac¬ 
teristic  length  (L),  corresponding  to  frequencies  greater 
than  a  crossover  frequency  a>co.  We  denote  the  Euclidean 
dimensionality  by  d,  and  call  the  density  dimensionality  of 
the  fractal  d  (the  so-called  Hausdor ff  dimensionality). 
(Thus  the  mass  increases  with  increasing  length  r  as 
r“).  The  exponent  giving  the  dependence  of  the  diffusion 
constant  on  distance  is  denoted  by  8,  such  that 

D(r)«r-*.  (1) 

Alexander  and  Orbach2  show  that  the  density  of  vibra¬ 
tional  states  in  the  regime  of  fractal  behavior  can  be  writ¬ 
ten  as 

Niaiocaj*-'  (2) 

where  d  is  the  fracton  dimensionality, 

j=2d/l2  +  d).  (3) 

In  Euclidean  space,  .\'(&>)a a>d~'.  The  crossover  frequen¬ 
cy  scales  as 


relevant  to  the  predictions  of  Ref.  2.  We  shall  describe 
the  relevance  of  these  experiments  to  fractal  behavior 
below,  but  first  it  is  necessary  to  make  explicit  the  cross¬ 
over  from  phonon  (long-length  scales)  to  fracton  (short- 
length  scales)  density  of  states.  In  particular,  we  need  to 
derive  the  normalization  coefficients  which  make  the 
above  Eqs.  ( 1 ) — (4)  quantitative. 

II.  PHONON  AND  FRACTON  DENSITY 
OF  STATES 

We  wish  to  normalize  the  long-wavelength  phonon  den¬ 
sity  of  states  to  a  volume  Ld.  This  is  because  the  phonon 
character  of  the  elementary  vibrational  excitations  ter¬ 
minates  at  the  (minimum)  length  scale  L.  Thus  we  set 

\p^co)  =  d(L/ald[iw^-'/(a)D}d]  ,  (5) 

where  a  is  an  atomic  distance  which  sets  the  shortest 
length  scale  in  the  problem  (i.e.,  fracton  behavior  is  ob¬ 
tained  for  length  scales  on  a  descending  basis  between  L 
and  a).  Thus 

tom  =  ia /L)u>d  ,  (6) 

where  ojd  is  the  apparent  Debye  frequency  as  projected  by 
the  low-frequency  velocity  of  sound.  The  integral  of  Eq. 
(5)  from  0  to  aico  equals  unity,  indicating  that  there  is  one 
mode  per  volume  Ld  as  required.  Thus  at  the  crossover 
frequency  one  finds 


<ocoocI- 2+9  /2  .  (4) 

Recent  experiments  of  Kelham  and  Rosenberg'  suggest 
that  experiments  involving  the  heat  capacity  and  thermal 
transport  of  epoxy  resins  may  be  exhibiting  behavior 


•Vph,wco)  =  <f/a>co  •  (7) 

The  fracton  regime  is  a  little  more  complex.  From  the 
definition  of  the  Hausdorff  dimensionality  there  will  be 
iL/a)*  atoms  per  molecule.  This  leads  to  an  appropriate 
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fracton  density  of  states.’ 

A'fr(w)=l(L/j)i[(cl>^-1/(£JFD|J]  •  (8) 

We  have  introduced  another  quantity,  the  so-called  frac¬ 

ton  Debye  frequency, 

ojpiL  /a)d/~  .  (9) 

This  quantity  sets  the  upper  limit  on  the  frequency  regime 
of  fracton  behavior.  Indeed,  if  we  integrate  Eq.  (8)  from 
<uc0  to  6>fd  we  find  I L  /a  )3—  1  modes  as  required.  That  is, 
there  is  one-phonon  mode  per  volume  Ld,  and  (L/a  )3—  1 
fracton  modes,  so  that  in  total  one  finds  ( L/a )3  modes  per 
volume  ad.  At  the  crossover  frequency,  Eq.  (8)  reduces  to 

.\(,i(oQQ)=d/ojco  .  (10) 

We  therefore  find  the  important  relationship  at  crossover, 

Arfr(<uco)/Arph(Wco)*j/d  .  (11) 

Because  we  believe  d  <  d,  we  shall  discover  that  the  exper¬ 

iment  appears  to  be  in  conflict  with  the  ratio  (11),  even 
though  the  remainder  of  the  spectrum  appears  to  be  con¬ 
sistent  with  the  ideas  of  fractal  behavior  (see  Note  added 
in  proof). 

III.  APPLICATION  TO  PERCOLATION 

Though  the  experimental  portions  of  this  paper  are  cer¬ 
tainly  not  described  by  percolating  networks,  it  is  of  in¬ 
terest  (and  possibly  of  experimental  importance)  to  carry 
through  the  ideas  of  the  previous  section  to  illustrate  an 
example  of  fractal  structure.  It  is  necessary,  however,  to 
distinguish  between  the  infinite  cluster  and  finite  clusters 
when  calculating  the  full  density  of  vibrational  states  for 
percolating  networks. 

As  shown  in  Ref.  2,  one  can  use  the  concept  of  fractals 
for  percolating  structures,  with  (L/a r  atoms  (sites)  per 
molecule,  where  L  is  the  percolation  correlation  length 
Denoting  the  infinite  cluster  by  the  superscripts  IC,  we 
find  Eq.  (8)  becomes 

A'}f(w)  = S(L/a)3[((o)3-'/(coFD)3]  .  (12! 

The  finite  clusters  need  to  be  considered  because  the 
(specific-heat)  density  of  states  will  be  the  sum  of  the  in¬ 
finite  and  finite  cluster  density  of  states.  The  probability 
per  site,  or  atom,  of  belonging  to  a  cluster  of  size  R  <  L  is 

[P(R)]Fc  =  (d-d)(R/a)''1+‘'W'  ,  (13) 

where  the  subscript  FC  means  finite  cluster.  Equation 
(13)  leads  to 

fL[P(R)]f:CdR=l-(a/L)'i-11 ,  (14! 

*  tf 

the  probability  of  belonging  to  a  finite  cluster.  By  defin¬ 
ing 

R'aji/o  d  ,  <  1 5 ; 

only  dusters  with  R  >Ruj 1  contain  modes  of  frequency 
cj.  Next,  the  total  number  of  sites  >or  atoms'  in  a  volume 
Ld  is,  from  Eq.  (14), 


(I  /a )dl  1  -ta/L  >d-‘ J  )FC-r  [\L/a>3)lc=*iL  /a )d  , 

(16) 

and  therefore  not  critical. 

Equations  (13)  and  (14)  together  enable  us  to  calculate 
the  total  number  of  modes  in  the  finite  cluster.  We  have 

A'fFrc(a>)  =  (L/fl)rf  fL  id/ojHcj/curoiipiRidR 

=  d  /(O 

X\[a/R(ci»)'3-d  -{a/L)2-d'\  .  (17) 

Adding  N}f(u)  from  Eq.  (12)  to  Eq.  (17),  and  using  Eq. 
(15),  we  find  the  remarkable  result  that  the  total  fracton 
density  of  states  on  a  percolating  network  equals 

Nu(o)=3(L/a)d(\/<1>)(u/o,fo)lx‘/,2+e'  .  (18) 

Should  we  integrate  Eq.  (18)  over  the  fracton  frequencies, 
we  would  find, 

Hlpjj 

f  Nh(o,)da)  =  (d/d)(Ld-\) ,  (19) 

or  only  d/d <  1  modes  per  atom,  and  not  one  full  mode. 
This  reflects  the  missing  center  of  mass  modes  of  the  fi¬ 
nite  clusters.  The  proof  follows  from  the  integration  of 
Eq.  (18)  from  qj(R),  the  smallest  frequency  allowed  for  a 
cluster  size  R,  to  the  maximum  fracton  frequency  <oD, 

CZNh{co)db,=R^ 1 

=R3[1-(1/R3)] 

=R3- 1,  (20) 

which  proves  the  statement  for  the  normalization  we  have 
used.  It  is  intriguing  to  speculate  that  this  "missing 
mode"  might  be  the  analogous  quantity  for  percolation 
systems  that  the  two-level  systems  are  for  amorphous  sys¬ 
tems. 

The  total  number  of  missing  modes  is  1  —  (d/d)  per 
atom  [see  Eq.  (19)].  Their  frequencies  are  unknown,  but 
their  mass  distribution  can  be  calculated.  There  are 
R  1  +d'  clusters  between  R  and  R  4  dR.  leading  to  a  mass 
distribution  of  mass  M~r  between  M(=R3)  and  M  +dM, 
with  r  =  (d/<7)  +  1.  The  missing  mode  frequencies  ft*  (tl 
for  “two-level"  systems,  and  M  for  the  mass  of  the  clus¬ 
ter)  must  be  less  than  [from  Eq.  (15)] 

ft X <R.v  2+91/2 «  M -[/3  .  (21) 

Summarizing  the  results  of  this  section,  Eq.  (18)  exhib¬ 
its  the  fracton  density  of  states  for  a  vibrational  network 
on  a  percolating  structure.  The  fracton  density  of  states 
at  crossover  remains  that  of  the  infinite  cluster: 

.Vfria>t(,)  =  .VjSwc„i  =  d/(yco  ,  (22 

but  the  slope  for  higher  frequencies  ' i.e..  within  the  frac¬ 
ton  regime  r-  proportional  to  cj:d  :~e  instead  of 
if3'  •-*  for  t!i?  infinite  cluster  alone.  This  is  because  of 
the  contribution  of  finite  clusters  to  the  fracton  density  of 
states.  In  addition,  the  integral  over  the  finite-cluster  den- 
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sity  of  states  shows  that  one  mode  per  cluster  is  missing. 
This  mode  can  be  attributed  to  the  center  of  mass  motion 
of  the  cluster,  and  may  be  analogous  to  the  "two-lei  e!  sys¬ 
tems"  attributed  to  amorphous  systems. 

We  also  note  in  closing  that  percolating  structures  can 
easily  be  achiesed  in  magnetic  systems  by  simple  dilution. 
A  previous  paper4  showed  that  the  diffusion  equation  not 
only  maps  on  to  the  vibrational  problem,  but  also  on  to 
the  linearized  spin-wave  problem  for  ferromagnetic  sys¬ 
tems.  Consequently,  for  randomly  diluted  ferromagnetic 
systems,  the  spin  excitations  would  cross  over  from  spin- 
wave-like  at  low  energies  to  spin-fracton-like  at  higher  en¬ 
ergies,  with  a  density  of  states  proportional  to  arJ _l 
for  the  infinite  cluster  alone,  and  to  -l  for  the 

sum  of  the  infinite  and  finite  clusters.  Either  specific  heat 
or  neutron-diffraction  studies  on  randomly  diluted  fer- 
romagnets  would  be  interesting  to  compare  with  these 
forms. 


IV.  RELEVANCE  OF  FRACTON  THEORY 
TO  EXPERIMENT 

We  have  already  suggested  that  recent  experiments  on 
the  specific  heat  and  thermal  conductivity  of  epoxy  resins 
by  Kelham  and  Rosenberg''  may  have  exhibited  fracton 
properties.  They  have  shown  (from  an  analysis  of  their 
specific-heat  measurements!  that  .V (oil  a.  o)Z  for  Hu>/kB  <  8 
K,  but  is  proportional  to  to  for  8  K  <fuo/kB  <  50  K. 
Their  figure  for  .Vital  appears  to  exhibit  a  discontinuity 
between  these  two  regimes,  though  the  analytic  form  they 
have  chosen  does  not.  This  change  of  slope  is  exactly 
what  one  would  expect  if  the  epoxy  molecules  were  exhib¬ 
iting  fractal  behavior.  The  "crossover  frequency”  is 
chosen  to  be  8  K.  and  corresponds  to  a  length  scale  of  30 
A.  about  the  length  of  the  epoxy  molecule  (diglycidyl  eth¬ 
er  of  bisphenol  .4).  For  stoichiometric  hardening  it  is  also 
the  distance  between  cross  links  of  hardeners.  Use  of  the 
same  number  of  hardener  molecules,  but  of  differing 
lengths,  did  not  change  the  crossover  ftequency  [Nicholls 
and  Rosenberg  (unpublished)]  and  the  crossover  length 
scale  remained  the  distance  between  cross  links.  This  is  to 
be  expected  if  the  epoxy  molecules  alone  are  exhibiting 
fractal  behavior.  Preliminary  evidence5  suggests  that  in¬ 
creasing  the  amount  of  hardener  (i.e.,  reducing  the  dis¬ 
tance  between  cross-links)  tends  to  raise  the  crossover  fre¬ 
quency.  If  one  associates  the  length  L  in  Eq.  (4)  with  the 
distance  between  cross  links,  then  this  effect  is  in  the 
correct  direction.  It  suggests  that  the  effect  of  the  cross 
links  is  to  restore  the  true  Euclidean  three-dimensional 
character  of  the  lattice  vibrations,  and  that  the  fractal 
behavior  is  to  be  associated  with  the  behavior  of  the  epoxy 
molecules  between  the  connections  with  the  hardener. 

The  observed  power  (linear)  for  the  vibrational  density 
of  states  above  the  hypothesized  crossover  frequency  sug¬ 
gests  d  —  2.  Unfortunately,  at  the  present  time  we  have  no 
independent  estimates  for  d  and  9  for  epoxy  resin.  We 
note,  however,  that  the  length  scale  can  be  changed  for  the 
epoxy,  either  by  changing  the  amount  of  hardener,  or  by 
using  different  epoxy  molecules,  which  allows  the  use  of 
Eq.  (4)  to  determine  6.  Further,  d  can  be  independently 


determined  from  x-ray  scattering  'see  the  method  of 
Stapleton  ei  alf  >.  It  would  be  of  great  interest  to  see  if  d 
(now  overdetermined1  is  consistent  with  such  measure¬ 
ments. 

There  remains  a  problem  between  the  analysis  of  the 
specific-heat  data  of  Kelham  and  Rosenberg'  and  the 
theory  developed  above  [specificajly.  with  Eq.  (11)].  Be¬ 
cause  d  <d,  and  with  6  positive,  d  <d.  so  that  the  fracton 
density  of  states  at  crossover  will  always  be  less  than  the 
phonon  density  of  states  at  crossover.  Thus  the  density  of 
states  can  only  suffer  a  drop  at  crossover,  whereas  the  ex¬ 
periments  appear  to  exhibit  a  rapid  rise  at  the  frequency 
which  we  have  interpreted  to  be  crossover.  We  are  at  a 
loss  to  understand  this  difference  in  behavior.  It  is  possi¬ 
ble  that  there  is  an  additional  feature  to  the  vibrational 
spectrum  of  epoxy  resins  in  the  fracton  regime  which  adds 
a  constant  to  the  density  of  states,  but  we  are  unaware  of 
its  origin.  As  can  be  seen  below,  the  analysis  of  the 
thermal  conductivity  is  also  consistent  with  fractal 
behavior  above  the  frequency  we  have  identified  as  cross¬ 
over,  so  there  is  some  compelling  character  to  the  evidence 
for  fracton  excitations  above  about  8  K  in  epoxy  resins. 
However,  the  inconsistency  with  the  direction  of  the 
discontinuity  remains,  and  makes  our  hypothesis  some¬ 
what  unsettled  (see,  however,  the  S'cie  added  in  proof  I. 

Perhaps  a  more  extraordinary  feature  of  the  measure¬ 
ments  of  Kelham  and  Rosenberg  is  their  finding  that  if 
one  analyzes  the  thermal  conductivity  in  the  usual  kinetic 
theory  manner,  the  mean  free  path  exhibits  a  precipitous 
drop  to  less  than  an  atomic  spacing  above  the  frequency 
corresponding  to  8  K  (or  the  length  of  30  A).  Such 
behavior  for  the  effective  mean  free  path  was  first  noted 
for  glasses  by  Zaitlin  and  Anderson.  They  used  a  nearly 
vanishing  mean  free  path  [region  C  of  their  Fig.  (5)]  as  a 
method  for  analyzing  the  “plateau  region"  found  for  the 
temperature  dependences  of  the  thermal  conductivity  of 
almost  all  glasses. 

A  mean  free  path  of  less  than  a  lattice  constant  certain¬ 
ly  suggests  localized  states  according  to  the  Ioffe-Regel 
rule.8  We  carry  the  fracton  picture  further  and  note  that 
Domany  el  al.  have  recently  shown’  that  essentially  all 
states  are  localized  on  a  bierpinski  gasket  (a  well-known 
fractal  geometry1).  In  general,  for  fractal  structures, 

Rammal  and  Toulouse10  have  recently  shown  that  the 
conductance  g(L)  scales  with  length  for  fractals  as 

giDccL3-1-9 .  (23) 

With  the  use  of  the  scaling  relation  of  Abrahams  et  al.,u 
their  quantity  0(g)  becomes 

p[g)=^~=d-2-9 ,  (24) 

al  I’LL 

or,  in  terms  of  the  fracton  dimensionality  d,  using  Eq.  (3), 

f)ig)=d[  1  —12/Ji]  .  (25) 

The  argument  of  Abrahams  et  al. 11  leads  to  localized 

states  for  0<O.  The  specific-heat  analysis  of  Kelham  and 
Rosenberg5  suggests  that  J— 1  =  1,  or  3=2.  Use  of  this 
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value  in  Eq.  (25)  results  in  /3(L)  = 0,  so  that  the  scaling 
factor  of  Abrahams  et  al.  is  consistent  with  localized 
behavior. 

For  localized  fracton  states  one  would  not  expect  any 
contribution  to  the  thermal  conductivity,  but  rather  only  a 
scattering  of  the  lower-energy  extended  phonon  states,  ex¬ 
actly  as  the  lower-energy  two-level  systems12  have  been 
shown  to  do  in  amorphous  systems.’  We  are  making  a 
distinction  here  between  the  two-level  systems  which  are 
also  found  in  epoxy  resin-’  (and  in  nearly  all  amorphous 
materials’1’)  and  the  fracton  modes.  The  density  of  states 
of  the  former  are  usually  added  to  the  extended  phonon 
density  of  states,  while  the  latter  replace  the  phonons  as 
the  fundamental  excitation  above  coco  (but  see  our  specula¬ 
tion  concerning  the  origin  of  the  "two-level  systems"  for 
percolating  networks).  Localized  fractons  would  then 
manifest  themselves  by  a  sharp  drop  (vanishing)  of  the 
contribution  to  the  thermal  conductivity  for  excitations 
with  <o>toco.  A  conventional  thermal-conductivity 
analysis  would  result  in  a  negligible  mean  free  path  in 
such  an  excitation  region.  The  property  of  the  epoxy  resin 
which  best  maintains  the  correctness  of  our  interpretation 
is  that  the  frequency  at  which  the  "mean  free  path”  be¬ 
comes  of  the  order  of  an  atomic  spacing  is  almost  precise¬ 
ly  the  frequency  at  which  the  density  of  states  suffers  a 
discontinuity — crossover  from  a  Debye-type  to  a  fracton- 
like  form. 

At  yet  higher  temperatures,  the  plateau  in  the  thermal 
conductivity  ends,  and  begins  to  increase  (e.g.,  above  ~  10 
K  for  epoxy  resin).  This  could  be  caused  by  classical  hop¬ 
pinglike  transitions  between  localized  fracton  states, 
analogous  to  Mott’s  variable  range  rate  hopping14  for  lo¬ 
calized  electronic  states.  We  have  not  as  yet  carried 
through  a  detailed  analysis  for  thermal  transport  in  such  a 
regime. 

The  fracton  hypothesis  suggests  that  materials  exhibit¬ 
ing  a  significant  temperature  spread  for  the  plateau  re¬ 
gime,  in  the  thermal  conductivity,  possess  a  characteristic 
length  considerably  larger  than  an  atomic  spacing  between 
which  vibrational  excitations  are  localized  (we  hy¬ 
pothesize,  fractonlike).  For  longer-length  scales,  one 
passes  into  the  normal  Debye-type  regime.  The  width  of 
the  plateau  region  will  increase  as  the  crossover  length  L 
increases. 

It  is  interesting  to  compare  the  results  and  analysis  of 
Kelham  and  Rosenberg3  on  epoxy  resin  with  that  of  Zait- 
lin  and  Anderson7  on  other  noncrystalline  materials  (a 
borosilicate  glass  and  a  polycarbonate).  As  recognized  by 
Kelham  and  Rosenberg,  Zaitlin  and  Anderson  were  the 
first  to  associate  the  plateau  in  the  thermal  conductivity 
of  almost  all  glasses  (see  Zeller  and  Pohl13)  with  a  sudden 
drop  in  the  phonon  mean  free  path  for  frequencies  above  a 
minimum  &>0  (the  crossover  frequency?).  This  is  clearly 
exhibited  in  their  Fig.  5,  region  C.  The  similarity  of  their 
analysis  for  noncrystalline  materials,  and  that  of  Kelham 
and  Rosenberg  for  epoxy  resin,  where  long-range  correla¬ 
tions  ( —  30  A)  are  expected,  is  striking. 

We  admit  that  it  is  a  little  puzzling  why  ordinary 
glasses  should  exhibit  such  a  length,  roughly  independent 
of  the  character  of  their  constituents.  However,  our 
analysis  would  not  be  the  first  to  argue  for  such  extended 


correlations.  For  example,  Morgan  and  Smith13  have 
made  similar  arguments  (even  to  lengths  up  to  1000  A), 
though  with  a  very  different  model  in  mind. 

It  is  interesting  to  note  that  as  pointed  out  by  Morgan 
and  Smith  the  plateau  temperature  width  is  much  larger 
for  amorphous  Se  than  for  most  glasses  iZeller  and 
Pohl13).  It  is  known  that  the  characteristic  correlation 
length  in  amorphous  Se  is  much  longer  than  in  most 
glasses.  We  have  already  noted  this  phenomenon  above 
for  epoxy  resin.  Finally,  Phillips16  has  also  argued  for  ex¬ 
tended  correlations  in  amorphous  systems. 

Recent  irradiation  experiments  on  crystalline  quartz 
also  tend  to  confirm  our  interpretation  of  fractal  behavior 
in  glasses.  Laermans  et  al. 17  have  shown  that  electron  ir¬ 
radiation  of  crystalline  quartz  does  not  produce  a  plateau 
region  in  the  thermal  conductivity.  However,  neutron  ir¬ 
radiation  does  generate  a  plateau  regime  with  a  width 
slightly  larger  for  larger  neutron  doses.18  X-ray  scattering 
by  Grasse  et  al.'9  shows  that  neutron  irradiation  causes 
amorphous  regions  of  diameter  20  A,  while  Grasse  et  al.20 
show  that  electron  damage  does  not,  hence  the  lack  of  a 
plateau  region  in  the  thermal  conductivity  for  the  latter. 
Finally,  Grasse  et  al.2'  point  out  that  the  size  of  the  amor¬ 
phous  regime  tends  to  grow  slightly  (20%)  with  increasing 
neutron  irradiation.  This  is  consistent  with  the  slight  in¬ 
crease  in  plateau  length  found  by  de  Goer  et  al.'*  as  a 
function  of  increasing  neutron  dosage. 

V.  SUMMARY  AND  CONCLUSIONS 

In  summary,  we  contend  that  the  thermal  properties  of 
epoxy  resin  and  glasses  can  be  understood  on  the  basis  of  a 
crossover  from  Debye-type  behavior  at  low  frequencies  to 
fracton  behavior  at  higher  frequencies.  The  density  of 
states  correspondingly  changes  the  exponent  of  its  power 
law  from  d  —  1  to  3 —  1  [Eq.  (2)].  At  this  same  crossover 
frequency,  the  vibrational  states  change  their  character 
from  extended  to  localized,  thereby  profoundly  affecting 
the  thermal  transport  and  serving  as  a  possible  explana¬ 
tion  for  the  extremely  small  mean  free  path  of  phonons  in 
this  energy  region  as  extracted  from  more  conventional 
analyses.  The  crossover  frequency  is  proportional  to  the 
inverse  of  the  length  scale,  according  to  Eq.  (4),  fracton 
behavior  expected  for  shorter  lengths.  Any  increase  in 
this  length  should  therefore  reduce  the  crossover  energy 
and  increase  the  width  of  the  plateau  region  measured  in 
thermal-conductivity  experiments. 

Note  added  in  proof.  The  form  we  have  used  for  the 
normalized  fracton  density  of  states  assumed  that  the 
force  constant  and  mass  scaled  smoothly  through  the 
crossover  regime.  Recent  work  of  P.  F.  Tua,  S.  J.  Putter- 
man,  and  R.  Orbach  [Phys.  Lett,  (in  press)]  suggests  an  al¬ 
ternative  picture:  a  drop  in  to  vs  inverse  length  scale  at 
the  crossover  length.  This  leads  to  a  jump  in  the  density 
of  states  (instead  of  a  drop),  going  from  the  phonon  to 
fracton  regimes.  Excellent  agreement  with  the  experimen¬ 
tal  results  of  Ref.  3  is  obtained.  Finally,  an  effective 
medium  approximation  for  .Y(o>  has  been  obtained  very 
recently  for  percolating  networks  by  B.  Derrida.  R.  Or¬ 
bach,  and  Kin-Wah  Yu  (unpublished).  It  strongly  sup¬ 
ports  the  assumptions  of  Tua  et  al.,  giving  additional  con- 
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fidence  to  a  fracton  interpretation  for  short-length  scale 
vibrational  excitations  in  amorphous  materials. 
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The  frequency  dependence  of  the  electrical  conductivity  is  examined  for  a  classical  hop¬ 
ping  model  of  a  random  one-dimensional  system  in  the  presence  of  a  superposed  static  elec¬ 
tric  field.  The  effect  of  the  field  is  taken  as  a  constant  bias  for  the  left-right  jump  rates.  A 
general  expression  is  derived  for  the  mean  velocity  and  frequency-dependent  conductivity. 
Explicit  evaluation  of  these  equations  is  given  for  correlated  and  uncorrelated  hopping  rates: 

( 1 )  in  general  for  high  frequencies  and  (2)  to  lowest  order  in  the  disorder  for  all  frequencies. 

In  the  latter  case,  an  initial  decrease  in  ai.a>)  for  very  small  a  is  found,  <7*  =a0— fl  ia)J, 
cr/sdioi1.  For  larger  frequencies,  the  conductivity  crosses  over  to  the  form 
cr(<i>):>sh0  +  b,ii<ii)l/2,  previously  calculated  by  Alexander  and  Orbach.  The  a„b,  are  con¬ 
stants  which  depend  on  the  strength  of  the  bias  and  the  randomness.  The  crossover  fre¬ 
quency  increases  with  the  bias.  In  addition,  the  variance  in  the  autocorrelation  function  is 
calculated  in  the  long-time  limit,  for  weak  disorder  in  the  symmetric  case.  It  is  shown  that 
fluctuations  do  not  significantly  affect  the  determination  of  this  quantity  under  these  condi¬ 
tions. 


I.  INTRODUCTION 

The  frequency  dependence  of  the  electrical  con¬ 
ductivity  in  one-dimensional  systems  with  symme¬ 
trical  random  hopping  rates  has  been  calculated  re¬ 
cently  in  a  series  of  papers.1-3  Recent  electrical 
conductivity  experiments  on  the  quasi-one- 
dimensional  conductor  quinolinium  di- 
tetracyanoquinodimethanide  (Qn(TCNQ)J  have 
been  analyzed  using  these  results  by  Alexander 
et  al  *  All  of  these  studies,  with  the  exception  of 
Ref.  3,  have  only  treated  the  case  of  symmetric  hop¬ 
ping  rates,  and  hence  were  only  applicable  in  the 
small  electric  field  regime. 

Recently,  attention  has  been  directed  to  the  ques¬ 
tion  of  transport  with  random,  but  biased,  hopping 
rates.3-7  These  calculations  have  explored  the  time 
dependence  of  the  mean  displacement.  Our  purpose 
is  to  calculate  the  frequency  dependence  of  the  con¬ 
ductivity  in  the  nonlinear  electric  field  regime.  We 
shall  find  a  result  very  different  from  the  symmetric 
case  at  the  lowest  of  frequencies,  crossing  over  to  a 
response  which  is  the  same  as  found  for  symmetric 
hopping  at  higher  frequencies.  This  crossover  has 
been  anticipated  in  Ref.  3. 


Our  purpose  here  is  to  formulate  the  problem  in 
such  a  manner  that  the  frequency  dependence  of  the 
conductivity  can  be  calculated  directly.  Section  II 
describes  the  formal  calculation  of  the  ac  conduc¬ 
tivity  in  a  uniform  dc  field.  Section  III  exhibits  the 
high-frequency  limit  for  the  conductivity.  The  solu¬ 
tion  for  weak  disorder  is  developed  in  Sec.  IV,  with 
figures  exhibiting  the  real  and  imaginary  part  of  the 
frequency  response  for  different  biases.  In  both 
Secs.  Ill  and  IV  we  treat  the  case  of  correlated  and 
uncorrelated  hopping  rates.  In  Sec.  V,  we  derive  an 
expression  for  the  variance  of  the  autocorrelation 
function  in  the  weak-disorder  limit,  and  show  that 
fluctuations  are  unimportant  in  this,  and  the  long¬ 
time,  limit.  Finally,  we  summarize  our  findings  in 
Sec.  VI. 

II.  FORMAL  SOLUTION  FOR  THE 
ac  CONDUCTIVITY 
IN  A  UNIFORM  dc  FIELD 

We  treat  the  same  master  equation  for  hopping 
transfer  as  in  Ref.  1,  but  with  asymmetric  hopping 
rates.  Let  F„U)  be  the  probability  that  site  n  is  oc¬ 
cupied  at  time  t.  Then,  for  near-neighbor  hops  only, 
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(1) 


where  Y„  +  |  is  the  hopping  rate  to  the  left  (from 
n  +  1  to  n  i  and  X„  _  |  is  the  hopping  rate  to  the  right 
(from  n  —  1  to  n).  In  the  presence  of  an  applied 
static  electric  field  9o<  and  a  small  oscillating  field 
9  at  frequency  <u, 


<8P0+*W 

exp 

ktT 

(90+9e,-)el 

— 

exp 

r  *.r 

tr. 


x„. 


i  • 


(2) 


Once  Q„  and  are  known,  we  may  find  V  from 
Eqs.  (5)  and  (4).  From  Eq.  (6b),  2«  =  0  whence 


V  = 


2  |  [*„ -aXH  . ,( 1  -Ee“>')][QH  +RmEe‘ -]  | 

n _ _ _ 

2Q- 


(7) 


Writing  the  velocity  in  the  form  V ~v  +£e'",o(ti>), 
one  identifies 


II 


IB 


(8) 


where  we  have  taken  the  potential  drop  to  be  the 
same  along  all  segments  of  equal  length  /,  and  where 
Wm  is  the  symmetric  hop  rate  between  n  —  1  and  n 
(i.e.,  the  hop  rate  in  the  absence  of  an  applied  static 
field).  Expanding  in  9, 


and 


a(<o)— 


2  [(*,  -a*.  _ ,  )R,  +aX,_  ] 

H 

2C. 


(9) 


Y,=a(l-£e'*'Ur._l ,  (3) 


where  E=29el/kBT  and  a=expi—29gel /kBT)  is 
termed  the  bias.  The  velocity  of  the  particle  at  time 
t  is  given  by 


2  **.<*> 

d  * _ _ 

dt  2?.(l) 


m 


2  (Xn  -Y,  )P„(i) 

IB 


(4) 


Expanding  P„(f)  in  a  power  series  in  the  reduced 
field  E, 


P'(t)~Q'+R'Ee‘*+  ■  ■  ■  ,  (3) 


one  can  derive  expressions  for  Q,  and  £„  from  Eq. 
(1).  One  finds 


^(aC.+i-CJ“Ar,_1(ae),-C._I)=0,  (6a) 

itoR9  =tXK{aRH  +  \  Rh  )—Xm^^<iRH  — £»_  |) 

-oX'Q'^+aX'.tQ'  .  (6b) 


These  are  our  two  principal  relationships  (the  first  is 
merely  a  restatement  of  the  result  obtained  in  Refs. 
3—7).  We  can  evaluate  Eq.  (8)  immediately  for  our 
model.  Expressing  the  Qm  as 

G,=  (10) 

r-o 

Equation  (6a)  is  solved  immediately.  One  finds 

This  is  the  same  result  as  in  Ref.  7,  if  we  note 
a=exp(  —  29oel/kBT)  and  <  1  /X0  >  =exp(  9gtl/ 
kBT)(\/W),  where  WH  is  the  symmetric  hop  rate 
between  n  —  1  and  n  in  the  absence  of  an  applied 
electric  field.  Inserting  into  Eq.  (11), 

u=2<l/IF>-1sinh(«’0e//k,D,  (12) 

a  well-known  result,7  and  one  which  reproduces 
Rice  el  al .*  for  the  case  of  a  regular  lattice.  We  now 
go  on  to  explicit  evaluation  of  the  conductivity,  Eq. 
(9),  in  the  high-frequency  (Sec.  Ill)  and  weak- 
disorder  (Sec.  IV)  regimes. 


in.  HIGH-FREQUENCY  LIMIT 


We  examine  the  case  when  the  frequency  co  is  much  greater  than  any  of  the  hopping  rates  Xn,  Y„.  Using  the 
relationship  Eq.  (6b),  one  expands  in  inverse  powers  of  o>  to  find  R„  in  terms  of  Q„ : 

R»  «  G„  + .  +  TT^rl"  WO*  + 1  -  « : +  .0.  +2-<*XmXK  _  ,Q„  +aX*Q,  + , -a 2X>- iQ* 

tM  yicj) 


1(0 


+a1X'_]X'QH  +  i+aXH  _  _2  Q„  _ , -aXl_  ,Q„  )+0 
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One  now  uses  the  expression  for  Qn,  Eq.  (10),  in  Eq.  (13)  to  obtain  the  explicit  form  for  A,.  Knowing  Q„  and 
/?„,  at  ail  is  known  from  Eq.  (9).  We  exhibit  two  forms  for  a(a>):  the  X„  uncorrelated,  and  the  X„  correlated. 
The  former  becomes 

autK„iaj)=ci<X)+^-[2aHX>2-a(a  +D(X2)+ai]-a)(X)/{\/X)) 

16) 

+  — l-r[2aHX)i~-iaHl+a){X)(Xi)+a(l+a)HX})-a{l-a)<X1)/(\/X)]+ 

(ico  r 


(14) 


This  expansion  is,  in  principle,  valid  as  long  as  <u  » (  X* ) 1  /*.  This  is  an  interesting  form,  in  that  some  unusual 
averages  enter  for  asymmetric  transfer  rates  (a  <  1).  For  the  symmetric  case,  a=  1,  and  one  recovers  the  high- 
frequency  form  for  <7(<o)  given  in  Ref.  1  [their  Eq.  (7.8)].  It  is  interesting  to  note  that  the  inverse  moment  of 
the  Xn  enters  at  high  frequencies  in  the  asymmetric  case.  This  suggests  nonuniversal  behavior  for  distribution 
of  class  (c)  (in  the  notation  of  Ref.  1)  even  at  high  frequencies,  whereas  the  behavior  is  purely  universal  in  the 
high-frequency  regime  for  symmetric  transfer  rates. 

The  correlated  case  has  the  following  high-frequency  conductivity: 


2  Uo /X,)af 

r-' 


if  6t»(  X“),/".  Because  of  the  translational  invariance,  one  can  replace  averages  like  (X0X,/Xp)  by 
(X'X,  +  ,/X,+,)  for  any  n. 

In  the  symmetric  limit  (a=  1)  for  short-range  (sr)  correlations,  it  reduces  to 


<rcor. „<*>)=<*) Jr(>-u2>]+"  - ,  a»»<jr->1/*. 


(16) 


IV.  WEAK-DISOROER  LIMIT 

We  pass  to  the  limit  of  weak  disorder  in  order  to 
obtain  a  solution  of  Eq.  (6b).  To  this  order  the  ex¬ 
pression  for  the  frequency-dependent  conductivity 
developed  below  [e.g.,  Eq.  (23)]  is  exact  for  all  fre¬ 
quencies  (see  the  Summary,  Sec.  VI,  for  further  dis¬ 
cussion  of  this  question).  We  set 


where  we  shall  work  only  to  linear  order  in  e.  Using 
Eqs.  ( 10’)  and  (18)  in  Eq.  (9),  one  finds  to  lowest  or¬ 
der  in  the  disorder,  for  uncorrelated  X„, 

<rmco,l(o)=[U-a)/A] 

x|(lA42)(a/9_l-/30)<e2> 

+ a/(  1  —a) 


~  =A  +en,  (<S  )W1«A  .  (17) 

In  exactly  the  same  manner  as  before,  we  can  write 
the  exact  solution  for  Q„, 

G„=M/<l-a>+  .  (10) 

p-  o 

We  must  now  evaluate  RH  in  terms  of  the  Q„  from 
Eq.  (6b).  To  carry  this  through  (to  any  order  in  e, 
2 HT".  K  =0),  we  expand, 


Rn  ~  2  &pe*+p 

P 

•  98 

7"  2  2  Yp.^n  +p^n  +f  7*  ’  *  (18) 


+  [a/(  1  — a)i4J](e2) |  ,  (19) 

where  it  remains  necessary  to  calculate  the  0,  de¬ 
fined  in  Eq.  (18).  The  algebra  leading  to  a  solution 
is  tedious,  and  we  merely  exhibit  the  results  here. 
For  p  >0. 


*  *<«-l)  . 

(A.-1  )U-a) 

(A.— a)2 


( 1  — a)(  A—  1  HA.2— a)  ’ 

and  for p%  - 1, 

a  ( A.  —  1  )2  k,+, 


0P=- 


( A  —  a)(  1  — a)(A.2— a) 


(20a) 


(20b) 


where  A  is  the  largest  root  of 
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X2  —  X(1  +a  -t-icoA )+a  =0  ,  (21) 


fi=a/X,  and  a  <  1.  In  particular,  for  purposes  of 
Eq.  (19), 


a(3_  |  —0o= 


a  a  -(  1  —  2A. 
1  —a  A.2  — a 


(22) 


leading  to  our  principal  result, 


a  < e 2)  a(X-l)2 

a— ("»“  7+7r"Xr-~ 


(23) 


The  frequency  limits  of  Eq.  (23)  are  interesting,  and 
illuminate  the  essential  physics  of  the  transport  pro¬ 
cess.  The  behavior  of  the  roots  of  Eq.  (21)  are  essen¬ 
tial  to  our  analysis.  If  one  attempts  to  make  a  low- 
frequency  expansion  of  X,  one  is  faced  with  a  com¬ 
parison  of  ( 1  —  a)1  to  coA  ( 1  +a). 

2X  =  1  A~icoA 

+[(l-a)2+2icoA(l+a)-(coA)2]l/2  .  (24) 


Our  analysis  of  this  ratio  is  based  on  a  comparison 
of  the  “drift”  distance  with  the  “diffusive”  displace¬ 
ment,  exactly  as  in  Ref.  3.  We  shall  show  that  the 
very  small  co  (long-time)  limit  is  dominated  by  the 
drift  of  the  particle  under  the  influence  of  the  elec¬ 
tric  field,  the  larger  a  domain  by  diffusion.  The 
former  will  result  in  a  diminution  of  crm — (co)  with 
increasing  frequency,  the  latter  in  an  increase.  The 
former  is  a  consequence  of  electric  field  driven  drift 
causing  the  particle  to  encounter  rare  very  small 
transfer  rates  which  diminish  the  conductivity.  It 
would  not  have  experienced  these  transfers  if  its 
motion  were  purely  diffusive  on  the  same  time  scale 
of  \/co.  The  latter  regime  is  associated  with  shorter 
and  shorter  root-mean-square  distances  as  co  in¬ 
creases,  leading  to  fewer  encounters  with  small 
transfer  rates. 

The  crossover  condition  can  be  obtained  in  exactly 
the  same  manner  as  in  Ref.  3  by  calculation  of  the 
relevant  lengths.  The  drift  velocity  has  already  been 
derived  in  Eq.  (11).  Approximating  (weak  disorder) 
<  1  /X)  =A,  in  a  time  l/<u,  the  particle  experiences  a 
net  average  displacement 


lM,=*{l-a)/(oA  ■  (25a) 

Likewise,  on  the  same  time  scale,  the  particle’s 
root-mean-square  diffusive  distance  is 

lM={Aca)-xn  .  (25b) 

The  ratio  of  the  square  of  the  two  is 

/5nft//ifl  =  [(l-a)2/dV]/[l/da»]  , 

or,  more  simply  ( 1  —  a  )2  / A  co  which,  for  small  1  —  a, 
is  precisely  the  ratio  Eq.  (24).  Hence,  the  crossover 


i 


in  frequency  behavior  can  be  traced  to  the  same  ori¬ 
gin  as  for  nonlinear  electric  field  effects  as  discussed 
in  Ref.  3. 

For  the  first  regime  (lowest  co ),  one  finds 


ouncm(co)~aQ-aico2+a2ico2  +  •  •  •  , 

co  <  ( 1  — cr  )2/A  , 

while  for  the  latter  regime, 


(26) 


<TUKOr(co)~b0  +  b{(ico){/2  +  •  •  •  , 

(27) 

w  >  ( 1  —a  )2/A  . 

The  coefficients  are  positive,  and  depend  on  the 
strength  of  the  disorder  through  Eq.  (23).  The  limit¬ 
ing  behaviors  are  easily  seen  if  one  plots  the  full 
solution  to  Eq.  (23).  Taking  A-  1  for  simplicity,  we 
have  calculated  the  frequency-dependent  part  of 
<j(co),  (A.—  1)2/(X2  — a),  against  co,  for  a=0.8  and 
0.9.  Figure  1  exhibits  the  real  part  of  this  factor  for 
the  interesting  frequency  regime  with  a =0.8.  One 
sees  a  crossover  at  co  = 0.046  between  the  two- 
frequency  regimes.  The  low-frequency  drift- 
dominated  limiting  regime  for  the  frequency- 
dependent  part  of  a  (co)  is  exhibited  in  Fig.  2,  and 
for  the  imaginary  part  in  Fig.  3.  Finally,  the 
higher-frequency  diffusion-dominated  regime- 
limiting  behavior  is  exhibited  in  Fig.  4.  Similar 
plots  for  a =0.9  are  exhibited  in  Figs.  5—7.  One  sees 
that  crossover  takes  place  at  lower  frequency, 
co — 0.01 1  from  Fig.  5.  This  follows  closely  the  pre¬ 
diction  of  Eq.  (24)  (as  it  must). 

These  curves  show  that  the  limiting  low- 
frequency  dependence  is  an  immediate  and  strong 
function  of  the  nonlinear  response.  For  small  elee- 


FIG.  1.  Plot  of  the  real  part  of  the  frequency- 
dependent  part  of  <7ico ),  (A.—  1  )V(A.:— a  I  (defined  in  the 
text),  vs  co  for  the  bias  parameter  a =0.8,  and  background 
regular  inverse  hop  rate  A  =  1 .  Note  the  crossover 
behavior  near  cj= 0.046. 
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trie  fields,  the  drift  regime  is  limited  to  very  small 
frequencies,  increasing  in  scale  as  the  electric  field  is 
enhanced.  For  the  symmetric  case,  a  =  1,  we  recover 
precisely  the  expression  for  eo)  given  in  Refs. 
2  and  9,  to  lowest  order  in  eo. 

For  correlated  Xn,  the  argument  goes  through  as 
before,  with  the  following  results  for  the  symmetric 
case  (a=  II: 


ocar(eo)t=\/A  +  , 


<«o> 


X-l 


X+l 


„ X-l  1  *  ,  *  1 


(28) 


and  the  asymmetric  case  (a  <  1): 


•• 


7T_ 


FIG.  2.  Same  as  for  Fig.  1,  but  for  small  to  on  an  ex¬ 
panded  scale.  The  straight  dashed  line  has  slope  2,  show¬ 
ing  that  £7*(k>)  falls  off  as  to1  for  smallest  to. 


,  ,  a  <4>  a  (X  —  1  )2  .  1  *  ,  v 


<X-a)2  af  .  a(X-l)2  1  ,,  ,  . 

+  — - —  —  (1— ala' 


X2— a  \* 


X2— a  X' 


V.  FLUCTUATIONS 
IN  THE  AUTOCORRELATION 
FUNCTION 

One  of  the  more  intriguing  questions  which  exact 
treatments  of  random  one-dimensional  systems  have 
been  unable  to  treat  is  fluctuations  about  the  ensem¬ 
ble  average.  The  simplest  quantity,  the  autocorrela¬ 
tion  function,  is  an  example  of  how  difficult  such 
problems  are.  As  discussed  in  Ref.  1,  Sec.  X,  the 
variance  of  the  autocorrelation  function  Poit),  in¬ 
volves  the  calculation  of 


(29) 


I - - - 

<f*$(t)> f'+‘~  da>'P0la>')Pol<»-a>')')  , 

(30) 

where  P0(eo)  is  the  Laplace  transform  ( Jf )  of  P0(t), 
where 

?„(«)«  f~e-*P„(t)dt. 

Because  the  P,(ot)  are  correlated1  for  different  eo, 
the  averaging  process  in  Eq.  (30)  does  not  commute 
with  the  convolution,  and  one  cannot  simplify  fur- 


! 


•9  '  I- 


FIG.  3.  Same  as  for  Fig.  1,  but  for  the  imaginary  part 
of  <7<a),  <7 ,(<o),  for  small  eo  on  an  expanded  scale.  The 
straight  dashed  line  has  slope  3,  showing  that  <r( (to)  in¬ 
creases  as  to'  for  the  smallest  to. 


FIG.  4.  Same  as  for  Fig.  1,  but  for  the  real  (R)  and  im¬ 
aginary  If)  parts  of  (X- 1  )2/( A2 — a)  above  the  crossover 
frequency  on  an  expanded  scale.  The  straight  line  has 
slope  y,  showing  that  <7/(&>)  rapidly  reaches  this  limiting 
behavior,  but  Ogla)  does  so  more  slowly. 
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FIG.  5.  Plot  of  the  real  part  of  the  frequency- 
dependent  part  of  <r(<u),  (X- 1  )J/(X2— a )  (defined  in  the 
text),  vs  oi  for  the  bias  parameter  a  =0.9,  and  background 
regular  inverse  hop  rate  A  —  1.  Note  the  crossover 
behavior  near  a =0.01 1 . 
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FIG.  6.  Same  as  for  Fig.  5,  but  for  on  an  ex¬ 
panded  scale.  The  straight  dashed  lines  have  slopes  3  and 
y,  for  comparison  with  limiting  behaviors  of  <T/(a>)  for 
frequencies  below  and  above  crossover,  respectively. 


ther  that  expression. 

Using  the  approximation  of  weak  disorder,  we 
have  been  able  to  make  some  progress.  The  Laplace 
transform  of  Eq.  (1),  with  the  initial  condition  that 
the  particle  is  at  the  origin,  n= 0  at  time  t—  0,  is 

tuF.(<u)=8^0+jr.[aF.  +  1(m)-F,(<u)] 

+*„-!  (31) 


As  before,  we  expand  in  the  e. : 

F.(tu)=C.(<u)+  2  a«.rer 

+  2  2  f -+•...  ,  (32) 

where  Q.,  a^r,  and  /3„,p  .  are  all  functions  of  a>.  A 
little  algebra  shows  that  if  the  X„  are  not  correlated. 


<?„(«)/*,,(«’))  —  <F,(w))<FB(<u')>  =  (e2)  2  ,  (33) 

fm-rn 

the  cancel  when  taking  the  difference  (33).  We  need  only  find  the  a^f(u>)  to  solve  Eq.  (33).  These  are 
found  as  before  by  use  of  Eq.  (32)  in  Eq.  (31).  One  finds  the  relations 


wC«=6ii.o+(l/^)(aGi,  +  i-G«  +Qn-i~aQ«  ) 

and 


(34a) 


««v  * ( 1  ,A )taa»  +  \.r  +“»  - t.j  ~ ( 1  +  a  ] 

-(l/Ai)S^(aQ^l-Qm)-(l/A2)6l,.Up(Qt.l-aQm) .  (34b) 

If  X  and  fi  are  the  two  roots  of 

aXz-X(l+a  -Mw)+1=0  ,  (35) 

with  |  X  |  >  \n\,  then  the  Q„  are  given  by 

n*0-.  Q^Akn^'Ak-n),  0-.  Q^Ak'  +  'fiAk-n) . 


The  aM  are  likewise  given  by 
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P> 0,  a%.= 


-  n>p  + 1 

(X-M>2 

«S, 


and  for 


— 1>  an,r~ 


(|-<‘v^,->,  nsp+l. 


(X-#<)2 
_ *fL 


(X-m>2 

Using  these  expressions,  Eq.  (33)  becomes 
< F0(a>)F0< *>' ) >  -  < F0(« ) ><■ ?(>(<»' )  >  =  < C2 > 


(X-l)(l-*i)X\  n  <p  . 


1  )2(Xj—  1  )2/iiP2+X|Xj(  1  —  #i|)2(  I  — P2)2] 

(X,  — /*|)2(A.2 — M2)2^l^-2  — 


(36) 


where  Xh/i,  and  X2,M2  are  the  roots  of  Eq.  (35)  (remembering  |  X  |  >  |  p  | )  for  w  and  <u\  respectively.  For  the 
symmetric  case,  a  =  1  and  Eq.  (36)  “simplifies”  to  (#at  =  1/X!t  p2=  1/X2), 


(P0(6>)P0(V')>-<P0(<U))(P0(V'))=<C2) 


2X|X2 


(X,+  l)2(X2+l)2(XfX|-l)  ' 


(37) 


We  wish  the  inverse  Laplace  transform  of  this  ex¬ 
pression,  a  formidable  task  which  so  far  has  eluded 
us  for  general  However,  for  small 

X=*  1  +  {Aai)ln,  whence  Eq.  (37)  reduces  to 

<  PQ(ai  )P0(ai' > )  —  (  Po(& ) )  (  F0(<u' ) ) 


if2) 

16^/A  ( v'oj 4. y/a' ) 


(38) 


FIO.  7.  Same  as  for  Fig.  S,  but  for  the  real  part  of 
(X—  1  )2/(X2-a)  above  the  crossover  frequency  on  an  ex¬ 
panded  scale.  The  straight  dashed  line  has  slope  y  for 
comparison  with  the  predicted  limiting  behavior. 


The  inverse  Laplace  transform  of  (v'^+v'^T')-'  is 
known,  so  that  the  long-time  behavior  of  the  vari¬ 
ance  of  P0lt)  is  approximately  given  by 

<ri(-)>-<r. 

M (39) 

t—»  00  . 

Comparing  with  (P0(r)>  =  v'A7(2v/irtl/J),  we  find 


<PiU))-(P0lt))2 

<P0(f)>2 


(c2) 
16vl  A2 


1/2 


(40) 


This  is  a  highly  satisfying  result,  for  it  demonstrates 
that  the  relative  variance  of  P0(t)  falls  off  with  in¬ 
creasing  time  as  t~ 1/2  for  the  symmetric  case  with 
weak  disorder.  The  asymmetric  case  is  more  com¬ 
plicated,  but  it  is  at  least  formulated  [Eq.  (36)]  in 
terms  of  the  Laplace  transformed  quantities.  The 
full  problem,  for  arbitrary  disorder,  seems  out  of 
reach  at  present. 


VI.  SUMMARY 

We  have  considered  the  random  one-dimensional 
near-neighbor  hopping  transport  problem  for  asym¬ 
metric  hopping  rates  with  constant  bias.  This  would 
correspond  to  the  problem  of  measurement  of  the 
frequency-dependent  conductivity  in  the  presence  of 
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a  large  dc  electric  field.  We  have  derived  expres¬ 
sions  for  the  drift  velocity  and  aito)  and  have 
evaluated  them 

(a)  exactly  for  the  drift  velocity  contribution  to 
the  dc  conductivity, 

(b)  exactly  in  the  high-frequency  limit,  and 

(c)  exactly  for  weak  disorder. 

In  case  (c),  we  have  shown  that  a  crossover  frequen¬ 
cy  exists,  proportional  to  the  bias,  (i.e.,  strength  of 
the  electric  field)  below  which  the  real  part  of  the 
conductivity  diminishes  as  and  the  imaginary 
part  increases  as  o?.  Above  the  crossover  frequency, 
the  real  and  imaginary  parts  increase  as  to  .  We 
have  exhibited  numerous  graphs  which  show  the 
real  and  imaginary  parts  of  a(a)  over  the  applicable 
frequency  range. 

We  have  limited  ourselves  to  the  first  correction 
term  for  the  Rn  [Eq.  (18)]  when  we  work  to  lowest 
order  in  the  disorder  (case  (c)  above].  To  that  order, 
the  frequency  dependence  we  have  found  for  <r(co)  is 
exact.  One  can  pose  the  question  of  the  effect  on 
aio>)  of  introducing  higher-order  corrections  in  Eq. 
(18).  Such  a  task  is  formidable,  but  we  know  the 
answer  in  the  symmetric  case  (a=  1).  There,  the 
aiu)  can  be  expanded  as  a  power  series  in  <<*>. 
One  finds9  that  these  terms  carry  with  them  only 
higher  powers  of  <o1/2: 

rfoj)<za{0)+CH<e2))d>>/2+0«e3))(o+---  . 

Thus,  to  lowest  order  in  (t2),  one  obtains  the  first 
frequency-dependent  correction  to  er(<y),  proportion¬ 
al  to  g>i/2.  To  higher  order  in  e,  one  obtains  yet 


higher  powers  of  uW2.  This  suggests  (but  we  have 
no  proof)  that  the  next  term  in  an  expansion  in  the 
disorder  for  the  asymmetric  case  would  yield  fre¬ 
quency  corrections  of  a  higher  power  than  we  have 
found  in  Sec.  IV.  Said  another  way,  it  may  be  the 
case  that  the  weak-disorder  limit  may  in  fact  gen¬ 
erate  only  the  lowest-frequency  correction  to  the  dc 
conductivity  in  the  asymmetric  limit  (as  it  does  in 
the  symmetric  case),  with  higher-frequency  correc¬ 
tions  arising  from  higher-order  terms  in  the  expan¬ 
sion  in  ,'owers  of  the  disorder  (as  it  does  in  the  sym¬ 
metric  case).  This  conjecture  remains  to  be  proven, 
but  if  true  it  means  that  our  expressions  for  case  (c) 
above  can  be  regarded  as  a  low-frequency  expansion 
for  the  situation  where  the  disorder  is  not  necessari¬ 
ly  weak. 

Finally,  we  have  shown  for  weak  disorder  that 
one  can  obtain  a  closed  form  expression  for  the  La¬ 
place  transform  of  the  variance  of  the  autocorrela¬ 
tion  function.  We  have  succeeded  in  obtaining  the 
inverse  for  the  symmetric  case,  and  have  shown  that 
fluctuations  above  the  mean  value  fall  off  as  t~'n 
for  long  times. 

ACKNOWLEDGMENTS 


The  authors  wish  to  acknowledge  very  helpful 
discussions  with  Dr.  J.  Bernasconi  and  Dr.  W.  R. 
Schneider,  and  receipt  of  their  interesting  work  on 
transport  for  asymmetric  hopping  rates  prior  to 
publication.  The  research  of  one  of  us  (R.O.)  was 
supported  in  part  by  the  U.S.  National  Science 
Foundation  and  the  U.S.  Office  of  Naval  Research. 


'Permanent  address:  Department  of  Physics,  University 
of  California,  Los  Angeles,  California  90024. 

*S.  Alexander,  I.  Bernasconi,  W.  R.  Schneider,  and  R.  Or- 
bach,  Rev.  Mod.  Phys.  Jl,  175  (1981). 

JS.  Alexander  and  R.  Orbach,  Physics  107B.  675  (1981). 

}S.  Alexander,  J.  Bernasconi,  W.  R.  Schneider,  and  R.  Or¬ 
bach,  in  Physics  in  One  Dimension,  Voi.  23  of  Springer 
Series  in  Solid  Stale  Sciences,  edited  by  J.  Bernasconi 
and  T.  Schneider  (Springer,  Berlin,  1981),  p.  277. 

V>.  Alexander,  J.  Bernasconi,  W.  R.  Schneider,  R.  Biller, 
W.  G.  Clark,  G.  Gruner,  R.  Orbach,  and  A.  Zettl. 
Phys.  Rev.  B  ]A,  7474  (1981);  S.  Alexander,  J.  Ber¬ 
nasconi,  W.  R.  Schneider,  R.  Biller,  and  R.  Orbach, 


Mol.  Cryst.  Liq.  Cryst.  S5.  121  U982). 

5M.  J.  Stephen,  J.  Phys.  C  14,  L 1077  (1981). 

6B.  Derrida  and  Y.  Pomeau,  Phys.  Rev.  Lett.  4£,  627 
(1982). 

7J.  Bernasconi  and  W.  R.  Schneider,  J.  Phys.  A  L 5,  L729 
(1982). 

!M.  J.  Rice,  S.  Strassler,  and  W.  R.  Schneider,  in  One 
Dimensional  Conductors,  edited  by  H.  G.  Schuster 
(Springer,  Berlin,  1975),  pp.  307—  309. 

’Reference  2,  and  J.  Machta,  Phys.  Rev.  B  24,  5260 
(1981);  R.  Zwanzig,  J.  Stat.  Phys.  2fl,  127  (1982);  I. 
Webman  and  I.  Klafter,  Phys.  Rev.  B  2 Si,  5950  (1982). 


I 


PHYSICAL  REVIEW  B 


VOLUME  38.  NUMBER  3 


1  AUGUST  1983 


Superconducting  transition  temperature  in  anodized  aluminum 
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We  have  measured  the  superconducting  transition  temperature  of  anodized  aluminum  films  of  gram 
sizes  ranging  from  less  than  100  to  3000  A.  The  transition  temperature  is  1  8  K  for  films  of  grain  size  100 
A  and  decreases  monotonically  with  increasing  grain  size  to  1  2  K  for  3000-A  grains  The  effect  depends 
only  on  the  volume  of  the  grains. 


The  enhancement  of  the  superconducting  transition  tem¬ 
perature  Tc  in  granular  aluminum  has  been  the  subject  of 
numerous  experimental  investigations. 1-3  An  unavoidable 
difficulty  with  the  granular  aluminum  is  (hat  (he  grain  size 
is  strongly  dependent  on  the  partial  oxygen  pressure  during 
evaporation,  making  it  virtually  impossible  to  distinguish 
between  the  effect  of  the  reduced  grain  size  and  that  of  the 
presence  of  the  oxide.  This  led  us  to  study  anodized  alumi¬ 
num:  both  effects  are  still  present,  but  they  are  much  more 
easily  separated.  We  find  that  in  our  anodized  aluminum 
films  of  grain  sizes  from  —100  to  —3000  A,  the  Tc 
enhancement  can  be  attributed  to  a  variation  in  the  grain 
size,  in  agreement  with  the  conclusions  reached  by  Chui 
et  at?  for  granular  aluminum  films.  Our  results  indicate 
that  the  effect  is  independent  of  the  surface-to-volume  ratio 
of  the  grains:  the  only  relevant  parameter  seems  to  be  the 
grain  volume.  All  samples  reported  here  were  fabricated  by 
evaporation  in  a  high  vacuum  and  further  processed  with 
the  use  of  photolithographic  and  anodization  techniques. 
The  substrates  were  precleaned  with  the  use  of  acetone 
washes  in  an  ultrasonic  cleaner  followed  by  an  alcohol  vapor 


FIG  l  Variation  of  Tr  with  sample  resistance  O.  O  series:  +, 
R  senes:  •.  M  series.  A.  P  series. 


bath.  Aluminum  was  evaporated  from  an  electron  beam 
gun  in  a  high  vacuum  (base  pressure  <3xl0-6  torn  onto 
heated  sapphire  substrates.  Substrate  temperature  varied 
from  40  to  250  °C,  resulting  in  a  variation  in  grain  size  from 
300  to  4000  A.4  The  film  thickness  was  measured  with  a 
crystal  quartz  monitor.  With  use  of  standard  photolitogra- 
phy.  the  films  were  etched  to  leave  samples  that  were 
—  0.28x1.22  mm.  The  thicknesses  of  the  samples  were 
further  varied  with  an  anodization  method.-1  A  constant 
current  was  passed  through  a  solution  of  saturated  boric 
acid  and  ammonium  hydroxide  with  the  use  of  a  platinum 
cathode  and  the  aluminum  sample  as  the  anode.  This  con¬ 
verted  the  aluminum  into  oxide,  resulting  initially  in  a 
linear  increase  in  the  anode-cathode  voltage  with  time  The 
anode-cathode  voltage  being  essentially  the  voltage  across 
the  oxide,  this  region  represents  an  increase  of  the  oxide 
thickness,  i.e..  a  thinning  of  the  aluminum  film  at  a  rate  of 
approximately  10  A/V.  At  a  film  resistance  of  50 R0.  where 
Ro  is  the  initial  film  resistance,  the  anodization  voltage  be¬ 
gins  to  increase  at  a  faster  rale,  indicating  that  the  film  is 
not  only  getting  thinner,  but  grains  are  being  isolated  from 
one  another.  In  this  manner  a  series  of  samples,  with  resis¬ 
tances  from  2  to  104  fi.  was  produced.  The  resistance  of 
these  samples  was  measured  in  a  liquid  4He  cryostat  by 
passing  a  constant  current  (7  —  10  ^A)  through  them  and 
measuring  the  voltage  drop  by  use  of  an  HP3455A  digital 


TABLE  I.  Properties  of  the  unanodized  samples. 


Sample 

series 

T, 

(C) 

1° 

(A) 

A 

Rjoo'R*: 

o 

40 

300 

230 

1.3 

R 

115 

800 

500 

2 

M 

250 

4000 

1000 

5 

P 

250 

4000 

800 

4 
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FIG  2.  Transition  temperature  as  a  function  of  the  cube  root  of  FIG  3  Transition  temperature  as  a  function  of  the  grains' 
the  gram  volume  a  O  series:  +.  R  series;  •.  V  series.  surface-to-volume  ratio  a  O  series;  +.  R  senes;  •.  M  series. 


voltmeter.  The  temperature  was  monitored  with  a  Lake- 
shore  Cryogenics  germanium  resistor  The  transition  tem¬ 
perature  Tc  was  defined  to  be  the  temperature  at  which  the 
resistance  falls  to  one-half  its  "normal”  value  above  the 
transition. 

For  low  normal-state  resistance  (R  <100  fl)  the  transi¬ 
tion  width  was  small  (  8777V  <  0.02 1.  At  high  resistances, 
however,  this  width  increases  and  for  one  sample  (series  O ) 
the  transition  was  incomplete  down  to  1.17  K,  the  lowest 
temperature  measured  in  these  experiments.  The  variation 
in  Tc  with  sample  resistance  is  shown  in  Fig.  1.  Four  dif¬ 
ferent  series  of  samples  are  plotted,  the  parameters  of  which 
are  listed  in  Table  I. 

In  Table  I,  T,  is  the  substrate  temperature  during  evap¬ 
oration.  t0  is  the  thickness  of  the  evaporated  films,  and 
Ajoo/A<2  is  (he  residual  resistance  ratio  of  the  unanodized 
films.  The  values  of  d0,  the  grain  size  of  thick  films  grown 
under  the  same  conditions,  are  not  known  exactly.  The 
values  listed  are  inferred  from  the  relationship  between  sub¬ 
strate  temperature  and  grain  size,4  and  are  consistent  with 
grain  sizes  observed  in  electron  micrographs  of  similar 
films.  The  error  involved  is  probably  ~  30%. 

Figure  1  shows  that  Tc  initially  increases  with  increasing 
resistance,  in  qualitative  agreement  with  the  measurements 
of  Deutscher  et  at. 1  for  the  case  of  granular  aluminum.  At 
approximately  103  fi  (the  resistance  region  in  which  anodi¬ 
zation  isolates  the  grains)  the  transition  temperature  levels 
off  and  then  decreases  for  even  larger  resistance  values. 
The  dependence  of  7V  in  log< 1  for  three  series  of  samples  is 
plotted  in  Fig.  2.  The  grain  size  dis  defined  as  d  -  ( d$t)'n\ 
i.e..  </  is  a  measure  of  the  grain  volume.  Only  the  samples 
with  resistance  values  R  <  IOORq  are  plotted,  where  Re,  is 
the  resistance  of  the  unanodized  film  of  the  series.  Since 
for  these  samples  do  is  virtually  unaffected  by  the  anodiza¬ 


tion,  the  sample  thickness  rcan  be  derived  from  the  sample 
resistance  at  room  temperature,  with  the  corrections  neces¬ 
sary  for  increased  surface  scattering.6  The  values  of  do  used 
for  Fig.  2  are  200,  900,  and  3500  A  for  the  O,  R.  and  M 
series,  respectively.  These  values  have  been  chosen  to  im¬ 
prove  the  agreement  between  the  three  series  and  are  well 
within  the  uncertainty  of  the  values  listed  in  Table  I. 

In  the  small  grain  size  region  id  <  200  A).  Fig.  2  agrees, 
within  experimental  error,  with  the  measurements  on 
granular  aluminum  reported  in  Ref.  3  For  large  grain  sizes 
the  slope  of  the  Tc  vs  logd  curve  appears  to  increase  (be¬ 
come  less  negative)  and  reach  bulk  TV  value  at  d  —  3000  A. 

One  of  the  models  proposed  to  explain  the  increased  Tc 
of  granular  aluminum  films  is  based  on  an  increase  in  the 
low-frequency  density  of  states  of  the  phonon  spectrum 
brought  about  by  the  high  surface-to-volume  ratio  of  the 
grains7  (“phonon  softening").  This  would  require  that  the 
increase  in  T,  be  proportional  to  S/F,  the  surface-to-volume 
ratio  of  the  grains.  For  our  grains  one  would  expect 
S/V  -ik/do  +  2/ 1.  where  k.  of  order  unity,  depends  on  the 
shape  of  the  grains.  A  plot  of  TV  vs  log(S/F)  is  shown  in 
Fig.  3.  For  any  one  series  of  samples  Tc  increases  mono- 
tonically  with  S/V,  in  agreement  with  the  phonon  softening 
model.  However,  for  every  sample  series  the  slope  on  the 
7V  vs  log(  S/F)  graph  appears  to  be  different.  We  conclude 
that,  as  evidenced  by  Fig.  2.  the  7V  enhancement  in  ano¬ 
dized  aluminum  films  is  a  function  of  the  grain  volume. 
The  phonon  softening  model  appears  to  break  down  for  thin 
grains. 
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Time  Decay  of  the  Remanent  Magnetization  in  Spin-Glasses 
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(Received  13  Jul>  1983.  revised  manuscript  received  1  November  1983) 

The  time  decay  of  the  thermoremanent  magnetization  (  otrm)  has  been  measured  in  1.0% 
Cu:Mn  and  2.6%  Ag:Mn  spin-glasses.  It  is  shown  that  (ttrm  is  neither  an  algebraic  nor  a  log¬ 
arithmic  function  of  time,  but  it  is  found  that  <ttrm  can  be  characterized  by  a  “stretched”  ex¬ 
ponential:  otrm”  troexpl  -  C(ut )1-*/(  1  -  /i) ].  Similar  time  dependences  appear  in  the 
disorder-difTusion  theory  of  Grassberger  and  Procaccia  and  the  cooperative-relaxation  theory 
of  Ngai.  but  neither  theory  in  its  present  form  is  directly  applicable  to  spin-glasses. 

FACS  numbers:  75.30.Kz.  75.50  Kj 


When  measured  in  a  small  static  field,  the  tem¬ 
perature  dependence  of  the  magnetization  of  a 
spin-glass  changes  abruptly  at  the  glass  temperature 
(Tg):  Above  Tg  the  magnetization  obeys  the  Curie- 
Weiss  law,  attributable  to  weakly  interacting 
paramagnetic  spins,  whereas  below  Tg  the  magneti¬ 
zation  is  nearly  independent  of  temperature,  indica¬ 
tive  of  the  spin-glass  state.  The  time  dependence  of 
the  magnetization  also  changes  dramatically  in  the 
vicinity  of  Tt :  In  the  paramagnetic  region  the  en¬ 
tire  magnetization  responds  rapidly  to  a  change  in 
field,  but  in  the  spin-glass  region  some  of  the  mag¬ 
netization  responds  much  more  slowly.1  One  tech¬ 
nique  of  investigating  this  viscous  behavior  is  to  ap¬ 
ply  a  Held  (H)  when  the  sample  is  in  the  paramag¬ 
netic  region.  Held  cool  it  through  Tg,  then  remove 
H  and  measure  the  Held-cooled  remanence  or 
“thermoremanent  magnetization”  (<rTRM).J 

The  exact  form  of  the  time  dependence  of  <rTRM 
has  not  been  previously  established.  Many  investi¬ 
gators3  have  reported  a  logarithmic  decay: 
o,TRM“°'o[l_(l~n)log(r)I,  where  oo  and  n  are 
constants.  Such  a  decay  is  unbounded,  however, 
and  must  be  merely  an  approximation,  valid  over 
some  Finite  interval  of  time.  The  Sherrington- 
Kirkpatrick  mean-field  model4  has  been  successful 
in  describing  many  of  the  observed  properties  of 
spin-glasses.  Calculations  based  on  this  model3  sug¬ 
gest  that  the  magnetization  should  decay  algebrai¬ 
cally:  <rTRM-<Vr1"*.  This  suggestion  seems  plau¬ 
sible  since  for  n  <  1  it  approaches  the  correct 
equilibrium  at  long  times  (zero  magnetization  in 
zero  Held6),  but  no  data  have  yet  been  published7 
supporting  an  algebraic  decay  for  <rTRM. 

We  have  made  magnetization  measurements  in 
the  interval  from  0.2  to  1000  sec  after  removing  H 
and  conclude  that  <ttrm  has  neither  an  algebraic  nor 
a  logarithmic  time  dependence.  We  do,  however. 
Find  that  the  time  dependence  is  accurately  charac¬ 


terized  by  a  “stretched”  exponential  of  the  form 

0’trm“0'oC*pI~  CUrl'-VO-n)].  (1) 

Here  the  exponential  factor  (C)  and  relaxation  fre¬ 
quency  (to)  can  be  chosen  to  be  independent  of 
temperature  throughout  the  spin-glass  region, 
whereas  the  prefactor  ( <r0)  and  time-stretch  ex¬ 
ponent  (n)  are  temperature-dependent  constants. 

We  have  made  time-decay  measurements  on 
three  different  samples:  1.0%  Cu:Mn,  2.6% 
Ag:Mn,  and  2.6%  Ag:Mn  +  0.46%  Sb,  ail  of  which 
show  qualitatively  similar  behavior.  Here  we  will 
present  only  the  2.6%  Ag:Mn  +  0.46%  Sb  data. 
Measurements  were  made  on  a  suck  of  thin 
(  —  25  fim)  foils  with  a  toul  mass  of  0.223  g.  The 
glass  temperature  for  this  sample  ( Tg  -  9.30  K)  was 
determined  from  the  maximum  in  the  magnetiza¬ 
tion  in  a  sutic  Field  of  3  Oe.  We  used  a  SQUID 
magnetometer  to  measure  a-TRM  as  follows:  (1)  a 
magnetic  Held  (H  -30  Oe  or  H  —  15  Oe)  was  ap¬ 
plied  to  the  sample  when  in  the  paramagnetic  re¬ 
gion,  (2)  the  sample  was  Held  cooled  through  T,  to 
a  temperature  in  the  spin-glass  region,  (3)  H  was 
removed  and  the  time  dependence  of  the  magneti¬ 
zation  was  recorded  by  a  computer-based  dau  ac¬ 
quisition  system;  (4)  after  1000  sec  the  remaining 
remanence  was  measured  by  warming  the  sample 
through  Tt  to  esubtish  the  base  line. 

Figure  1(a)  is  a  plot  of  crTRM  vs  log( /)  at  four 
temperatures  within  the  spin-glass  region.  The  fact 
that  the  slope  changes  with  time  shows  that  <rTRM 
does  not  decay  logarithmically.  Figure  1(b)  is  a 
log-log  plot  of  the  same  dau,  here  demonstrating 
that  otRM  docs  not  decay  algebraically.  In  Fig.  2  we 
show  that  <7trm  decays  by  a  stretched  exponential 
function  of  time.  The  time-stretch  exponent  ( n)  is 
most  easily  determined  by  plotting  log[-(<// 
dt)(ln<rTRM))  as  a  function  of  log(r).  This  is  done 
in  Fig.  2(a)  where  a  Fit  to  the  data  after  5  sec  yields: 
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FIG.  1.  (a)  Semiiog  and  (b)  log-log  plot  of  che  ther- 
moremanenc  magnetization  (<tt*m)  of  2.6%  Ag.'Mn 
+  0.46%  Sb  as  a  function  of  time  at  T/T, -0.771,  0.856, 
0.897,  and  0.966,  from  top  to  bottom,  respectively.  The 
solid  curves  are  the  best  stretched-exponential  fits  to  the 
data. 


» -0.694.  0.740,  0.766,  and  0.831  for 

r/T,- 0.771,  0.856,  0.897,  and  0.966,  respectively. 
The  stretched  exponential  nature  of  o-trm  >s  veri¬ 
fied  in  Fig.  2(b)  by  the  linear  dependence  of 
1o8(o>>trm)  on  tl~".  The  quality  of  agreement  is 
also  shown  in  Fig.  1,  where  we  plot  the  best 
stretched  exponential  fit  to  each  set  of  data.  The 
deviation  before  5  sec  can  be  attributed  to  the  decay 
of  induced  eddy  currents  in  the  metallic  samples  as 
is  evidenced  by  a  measurable  magnetic  absorption,* 
even  in  the  paramagnetic  region,  at  frequencies 
above  0.2  Hz.  After  5  sec,  <rTRM  is  accurately 
characterized  with  the  four  parameters  o-o,  C,  w, 
and  n.  Note  that  no  adjustment  is  allowed  in  the 
base  line  of  the  data  since  it  is  unambiguously  es¬ 
tablished  in  step  4  of  our  procedure. 

The  temperature  dependence  of  the  prefactor  and 
the  time-stretch  exponent  is  shown  in  Fig.  3.  At 
low  temperatures  (r<0.75rf)  cr0  decreases 
linearly  with  increasing  temperature  whereas  «  =*  7 
independent  of  temperature.  For  T  >  0.75  7",,  n  in¬ 


FIG.  2.  Determination  of  the  time-stretch  exponent,  it. 
(a)  Log-log  plot  of  -rfln(aTRM)/<ft  as  a  function  of  time. 
The  slope  gives  -  n  and  the  r  —  1-sec  intercept  gives 
Ccv1"*.  (b)  Semilog  plot  of  <tt»m  as  a  function  of  t'~‘. 
The  solid  lines  are  the  best  fits  to  the  data. 


creases  while  <ro  decreases  more  rapidly  than  at 
lower  temperatures.  The  temperature  dependence 
of  n  allows  us  to  extract  the  relative  values  of  <■>  and 
C  by  plotting  log(Coi1_"),  obtained  from  Fig.  2(a), 
as  a  function  of  1  -  ft.  This  is  done  in  Fig.  4  where 
linear  regression  to  the  data  yields  C  —  0.59  ±  0.05 
and  cv-(3  ±  l)xl0~6  sec-1.  The  only  difference 
we  measure  between  cooling  the  sample  in  //- 15 
Oe  and  in  H  -  30  Oe  is  a  linear  dependence  o0  on 
H.  it,  C,  and  <u  do  not  depend  on  the  cooling  field 
for  these  relatively  small  values  of  H\  thus,  the 
changes  in  n  and  <ro  for  T  >  0.75  f,  cannot  be  attri¬ 
buted  to  the  saturation  of  <ttrm  H-9 
Two  separate  theories  of  relaxation,  which  have 
recently  appeared  in  the  literature,  give  a  stretched 
exponential  time  behavior  similar  to  that  of  Eq.  (1). 
Neither  theory  has  yet  been  applied  to  spin-glasses, 
so  we  can  only  point  out  the  possibility  of  a  connec¬ 
tion.  The  first  theory,  due  to  Grassberger  and  Pro- 
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FIG.  3.  (a)  The  temperature  dependence  of  n. 

(b)  The  temperature  dependence  of  the  prefactor  (era) 
The  solid  circles  are  for  a  30-Oe  cooling  Held  and  the 
open  circles  are  for  H- 15  Oe.  The  dashed  lines  are 
guides  for  the  eye. 


caccia,10  considers  the  diffusion  of  particles  through 
a  d-dimensional  space  interspersed  with  randomly 
distributed  traps.  A  connection  to  spin-glasses  can 
be  made  by  using  the  model  of  Bantilan  and  Pal¬ 
mer11  in  which  the  energy  of  a  spin-glass  is  pictured 
as  a  labyrinthine  function  in  spin-configuration 
space  containing  many  maxima  and  local  minima, 
and  several  quasidegenerate  ground-state  minima. 
The  configuration  of  a  group  of  spins  travels  ran¬ 
domly  through  configuration  space  until  it  is 
trapped  into  one  of  the  ground-state  minima. 
Grassberger  and  Procaccia  find  that  the  number  of 
untrapped  “configurations"  (Ne)  is  given  by 

Nt  =*f(t)txpl -g  («)(«/  -»)]. 

where  the  time-stretch  exponent  is  related  to  the 
dimension  of  the  diffusive  space  by 
n-l-rf/(d  +  2).  Although  this  interpretation  is 
interesting,  three  experimental  facts  emphasize  the 
need  for  further  development  before  this  model 
may  apply  to  spin-glasses.  First,  otrm  obeys  a  sim¬ 
ple  stretched  exponential  of  time;  we  do  not  mea¬ 
sure  any  time  dependence  to  the  prefactor  I/(t)J. 


FIG.  4.  Semilog  plot  of  Cw'~*  as  a  function  of  1  -  n 
for  H  —  30  Oe  (solid  circles)  and  H  —  15  Oe  (open  cir¬ 
cles).  The  solid  line  is  the  best  fit  to  the  data.  The  slope 
of  this  line  gives  u-  (3  ±  l)x  10”‘  sec'1  and  the  inter¬ 
cept  gives  C  -0.59  ±0.05. 


Second,  we  can  choose  <•>  so  that  the  exponential 
factor  If  («)  I  is  independent  of  n.  And  third,  we 
find  n  &  j,  which  implies  a  nonphysical  (d  « 1 ) 
configuration  space. 

The  second  theory  giving  a  stretched  exponential 
time  decay  is  due  to  Ngai.12  The  Ngai  theory  treats 
the  cooperative  relaxation  of  a  primary  system  of 
dipoles  perturbatively  coupled  (coupling  constant 
V)  to  a  secondary  continuum  of  low-energy  excita¬ 
tions  whose  density  of  available  levels  is  linear  in 
energy:  N(E)-aE.  In  the  Ngai  model  the  mea¬ 
sured  susceptibility  of  a  sample  is  due  entirely  to 
the  particular  microstate  of  the  dipolar  system,  but 
the  rate  at  which  the  dipoles  approach  equilibrium 
is  influenced  by  the  continuum.  If  the  coupling 
could  be  ignored,  the  dipoles  would  make  transi¬ 
tions  at  some  time-dependent  rate,  1/ 7~0.  But  if  this 
coupling  is  not  negligible,  the  dipole  transitions  will 
excite  the  continuum,  which  in  turn  will  influence 
the  transition  rate.  Ngai  gives  the  susceptibility  as  a 
function  of  time  from  which  the  magnetic 

response  to  the  removal  of  H  at  f -0  is  easily 
found: 


MU)-H  fait')*' 


-Hi|i  o«xp 


,i-» 


1  —  n 


(2) 
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where  Et  is  a  cutoff  energy  for  the  low-energy  exci¬ 
tations  and  the  time-stretch  exponent  is  given  by 
If  Ec  and  1/7*0  are  set  equal  to  w,  then 
Eq.  (2)  has  the  same  form  as  Eq.  (1).  The  Ngai 
theory  gives  a  specific  value  for  the  exponential  fac¬ 
tor  in  terms  of  Euler's  constant:  C  —e~y  — 
0.5615...  which,  to  within  experimental  accuracy, 
is  the  value  we  measure.  The  Ngai  model  empiri¬ 
cally  characterizes  the  observed  time  dependence  of 
<rTRM.  but  a  significant  inconsistency  still  exists:  Eq. 
(2)  is  an  approximation  supposedly  valid  only  at 
long  times  (£e/»l),  which  for  £f-a>  ~  10-5 
sec”’  is  never  achieved  in  our  measurements.  We 
cannot  yet  explain  why  Eq.  (2)  seems  to  be  valid 
throughout  the  time  regime  of  our  measurements. 
Furthermore,  the  source  of  the  low-energy  excita¬ 
tions,  why  1 1 T0  should  equal  Ec,  and  the  tempera¬ 
ture  dependences  of  <r0  and  n  are  not  yet  under¬ 
stood. 

In  conclusion,  we  have  shown  that  the  decay  of 
<ttrm  in  spin-glasses  is  neither  a  logarithmic  nor  an 
algebraic  function  of  time.  It  is,  however,  accurate¬ 
ly  characterized  at  all  temperatures  within  the  spin- 
glass  region  by  a  stretched  exponential  (Eq.  (1)1 
with  four  adjustable  parameters.  The  prefactor 
(<r0)  and  time-stretch  exponent  (n)  are  tempera¬ 
ture  dependent,  whereas  the  exponential  factor  ( C) 
and  relaxation  rate  (<u)  can  be  chosen  to  be  in¬ 
dependent  of  temperature  throughout  the  spin-glass 
region.  In  addition,  we  find  o0  to  depend  linearly 
on  the  cooling  field,  whereas  n.  C,  and  w  are  in¬ 
dependent  of  H  for  H  ^  30  Oe.  We  point  out  pos¬ 
sible  connections  to  two  recent  theories  of  relaxa¬ 
tion,  but  emphasize  that  further  development  is 
necessary  before  either  theory  may  be  applied  to 
the  decay  of  aTRM  in  spin-  glasses. 
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A  brief  list  of  current  areas  of  research  in  high  field  physics  is  presented  covering 
most  of  the  presentations  at  this  Symposium.  Mora  detailed  description  is  given  for 
three  topics  for  which  high  magnetic  fields  are  required,  and  which  posseaa  unusual 
incerest.  These  are:  1)  Density  of  states  for  vibrational  states  on  fractals 
"Fractons",  2)  Thermodynamic  properties  of  exchange  enhanced  systems,  and  3)  p-atate 
pairing  in  thin  film  or  layered  superconductors. 


1.  INTRODUCTION 

This  Symposium  follows  at  least  two  others  in 
the  rapid  development  of  high  magnetic  field 
physics.1’2  In  addition,  a  survey,  now  rather 
aging,  of  opportunities  for  research  in  high 
magnetic  fields  has  been  prepared.2  The  pur¬ 
pose  of  the  present  paper  is  to  briefly  classi¬ 
fy  the  character  of  those  papers  to  be  pre¬ 
sented  at  this  Symposium,  and  then  to  describe 
in  outline  form  three  areas  which  are  of  par¬ 
ticular  interest  to  the  author. 

The  general  areas  of  research  in  high  magnetic 
fields  to  be  discussed  at  this  Symposium  can 
very  roughly  be  titled  as: 

1)  Collective  phenomena  (e.g.,  p-wave  super¬ 
conductivity) 

2)  Magnetic  structures  (e.g.,  phase  transi¬ 
tions,  magnetic  saturation) 

3)  Atomic-like  states  (e.g.,  exclton  structure 
and  dynamics) 

4)  Diamagnetism  (e.g.,  orientational  ordering 
of  large  molecules) 

5)  Thermodynamic  properties  (e.g.,  field 
dependent  susceptibilities) 

6)  Transport  properties  (e.g.,  quantum  oscil¬ 
lations) 

7)  High  energy  density  of  states  (e.g.,  vibra¬ 
tions  on  a  fractal) 


2.  DENSITY  OF  VIBRATIONAL  STATES  ON  A  FRACTAL, 
"FRACTONS" 

Fractals  are  open,  self  similar  atructurea, 
with  interesting  properties  as  a  function  of  the 
length  scale. ^  A  specific  example  would  be  a 
percolation  arrangement  where  the  number  of 
sites  on  the  infinite  cluster  (p  >  pc>  where  pc 
is  the  percolation  threshold  concentration) 
increases  not  as 

d 

r 

where  r  is  the  distance  and  d  the  Euclidean  di¬ 
mensionality,  but  rather  as 


vhere  d  is  an  effective  dimensionality,  equal 
to  d  -  (B/v)  in  terms  of  the  usual  percolation 
exponents. 2  This  behavior  occurs  for  short 
length  scales  in  comparison  to  the  coherence 
length  for  percolation,  $  .  For  larger  lengths 
one  finds  usual  Euclidean*1  properties.  If  now 
one  examines  diffusion  along  the  Infinite  clus¬ 
ter,  the  "dead  ends"  cause  a  length  dependence 
for  the  diffusion  constant: 

D(r)  «  r_lS 

where  again  for  percolation  £  •  (t-B)/v,  t  being 
the  conductivity  exponent. 


No  list  is  complete,  but  this  can  serve  as  a 
rough  outline  of  topics  unique  to  high  magnetic 
field  research. 

This  paper  will  explore  three  of  the  seven  areas 
listed  above.  The  remaining  four  will  be  well 
covered  by  others  at  this  Symposium.  Only  one 
of  these  three  represents  original  work  by  this 
author.  However,  the  significance  of  the  other 
two  warrants  soma  attention. 

Each  of  the  three  topics  is  described  below  in 
terms  of  the  physical  ideas  which  have  been 
developed,  and  the  possible  experimental  probes. 
Space  limitations  require  that  the  reader  be 
referred  to  the  original  treatments  for  the 
complete  details. 


The  diffusion  problem  along  a  fractal  can  be 
solved,  leading  to  the  ensemble  averaged  auto¬ 
correlation  function6 


<P0(t» 


-d/ (2+6) 


(1) 


where  the  particle  has  assumed  to  have  been 
localized  at  che  origin  at  time  t  ■  0. 


One  now  notes  that  the  form  of  the  diffusion 
equation  (Master  Equation)  is  the  same  as,  for 
example,  the  harmonic  vibrational  problem,  with 
a  simple  replacement  of  the  first  time  deriva¬ 
tive  by  the  second.  This  mapping  allows  us  to 
regard  the  lnverae  Laplace  transform  of  Eq.  (1) 
as  the  lattice  vibrational  density  of  states 
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(with  a.2  replacing  tha  Laplace  transform  spec¬ 
tral  parameter  e)  for  a  fractal  arrangement  of 
masses  and  springs.  Ona  finds 

N(w)  «  J\  p  -  (2d/(2+S)  ]  -  1  (2) 

For  Euclidean  systems,  p  *  d-1,  so  we  era  led 
to  define,  for  mode  counting  purposes,  a  recip¬ 
rocal  space  of  effecciva  dimensionality 

3  -  2d/ ( 2+3 )  .  (3) 

Ue  refer  to  these  states,  when  quantized,  as 
"fractons."  Their  properties  are  most  interest¬ 
ing.  Before  we  outline  them  in  more  detail, 
some  experimental  examples  are  of  interest. 

Our  attention  to  this  problem  was  aroused  by  the 
work  of  Stapleton  et  al.7  who  measured  ehe  spin- 
lattice  relaxation  time  for  low-spin  Fe(3+)  in 
three  hemoproteins.  These  large  molecules  were 
shown  by  x-ray  measurements  (counting  the  in¬ 
crease  of  the  number  of  alpha  carbons  with  dis¬ 
tance  for  myoglobin  at  230  K)  to  yield  a  value 
for  d  -  1.67  ♦  0.04,  certainly  not  integral. 
Their  data  for  the  spin-lattice  relaxation  time 
as  a  function  of  temperature  for  myoglobin 
azide  (MbN^)  are  copied  below: 


2  3  4  6  10  15 
T(K) 


Fig.  1.  The  electron  spin  relaxation  rate  of 
low-spin  Fe(3+)  in  MbNj.  The  race  is  fitted  to 
the  sum  of  a  direction  process,  varying  as  T, 
and  a  Raman  process,  with  taaperaturs  exponent 
6.29. 


Their  Interpretation  relied  on  the  use  of  the 
usual  two-phonon  integral  for  the  Raman  process 
relaxation  rate,  the  integrand  being  propor¬ 
tional  to  the  square  of  the  vibrational  density 
of  states.  Keeping  all  other  factors  tha  same 
as  for  Euclidean  space,  they  extracted  the 


exponent 

p  -  0.65  i  0.04  .  (4) 

They  did  not  report  ocher  data  which  would 
enable  ua  to  obtain  an  independent  estimate  for 
4.  The  use  of  self  avoiding  random  walks  as  a 
model  for  these  proteins  is  inappropriate.  For 
such  systems  in  d  »  3,  ® 

d  -  5/3,  4  -  4/3,  (5) 


leading  to  d  *  1  (p  *  0),  representing  one 
dimensional  vibrational  behavior. 

The  essential  condition  for  application  of  these 
ideas  to  physical  systems  is  that  tha  length 
scale  be  leas  than  the  Euclidean  correlation 
length.  For  lattice  vibrations,  this  implies 
that  the  frequency  be  greater  than  a  crossover 

frequency,  u  ,  itself  related  to  4  by  Che 
c.o . 

following  expression6: 


2 

W 

c.o. 


>  l-(2+5) 

“V 


(6) 


where  L  is  the  size  of  the  fractal  object  (e.g., 
tha  percolation  correlation  length,  or  the  size 
of  the  molecule)  in  units  of  the  monosier  length, 
and  the  frequency  scale  is  that  of  the  Debye 
frequency  appropriate  to  the  frectal  object. 

For  example,  Stapleton  at  al.  state  that  the 
temperature  range  1  -  20  K  is  associated  with 
wavelengths  of  from  10  to  10 3  bonds.  For  a 
large  molecule,  this  would  certainly  be  consis¬ 
tent  with  the  requirement  for  fractal  behavior. 


Thera  are  other  properties  of  fractal  vibration¬ 
al  states.  For  example,  the  vibrational  eigen¬ 
functions  are  local  and  should  not  contribute  to 
the  thermal  conductivity.  This  behavior  (though 
not  with  an  identification  of  fractal  proper¬ 
ties)  has  recently  been  reported  by  Kelham  and 
Rosenberg  for  epoxy  resin,  for  the  energy  range 
of  8  -  50  K  (their  measurements  spenned  the 
range  of  0.1  -  80  K).9 


It  is  clear  chat  the  identification  with  fractal 
behavior  depends  on  the  condition  (6),  which 
than  leads  to  a  vibrational  density  of  states 
(2).  The  experimental  consequences  are  Immedi¬ 
ate.  The  one  phonon,  or  direct  relaxation 
process  rate,  is  directly  proportional  to  the 
vibrational  density  of  states.  If  one  performs 
an  electron  spin  lattice  relaxation  time  mea¬ 
surement  at  sufficiently  high  magnetic  fields, 
it  is  possible  that  one  can  obtain  a  direct 
measurement  of  the  fracton  density  of  states. 

The  field  must  be  sufficiently  high  chat  the 
condition  (6)  is  satisfied.  Then  the  field  de¬ 
pendence  of  the  relaxation  rate  will  give  the 
energy  dependence  of  the  density  of  fracton 
states,  and  hence  a  value  for  p  using  Eq.  (2). 
The  crossover  magnetic  field  will  separetaly 
give  an  estimate  of  4  using  Eq.  (6).  That  is, 

5  >e  not  a  free  parameter,  in  that  it  is 
determined  by  the  crossover  behavior.  Finally 
the  factor  d  can  also  be  determined  from  x-ray 
measurements,  over-determining  all  the  fractal 
parameters. 


On  the  Physics  of  High  Magnetic  Fields 


5 


There  are  ocher  interesting  consequences  of 
fractal  behavior.  The  eigenstates  are  supposed 
to  be  localized.  This  could  lead  to  rather 
interesting  magnetic  resonance  bottleneck  ef¬ 
fects  in  that  the  spatial  transfer  of  excita¬ 
tion  will  be  diffusive  rather  than  wavelike. 

This  might  lead  to  strong  bottleneck  conditions 
for  both  the  direct  and  resonance  relaxation 
processes.  Here  too,  large  magnetic  fields 
would  be  useful  to  unravel  Che  dynamical  proper¬ 
ties  of  a  bottleneck.  For  example,  if  a  bottle¬ 
neck  is  found  for  the  resonance  relaxation  pro¬ 
cess,  the  field  dependence  of  the  strength  of 
the  bottleneck  will  give  a  direct  estimate  of 
the  fracton  lifetime  (the  analysis  is  similar 
to  that  of  Gill,10  but  using  fractons  Instead 
of  phonons I, 

3.  THERMODYNAMIC  properties  of  exchange 
ENHANCED  SYSTEMS 

The  effect  of  magnetic  fields  upon  the  thermo¬ 
dynamic  properties  of  Fermion  systems  [e.g., 
the  nonlinear  magnetic  susceptibility  and  the 
field  dependent  specific  heat],  depends  on  the 
relation  between  field  H  and  Fermi  temperature 
Tp.  Even  for  extreme  fields,  this  ratio  is 

small  (100  Tesla  is  the  equivalent  of  170  K  in 
Zeeman  splitting).  Exchange  enhancement  caused 
by  electron-electron  interactions  in  a  matal 
can  significantly  enhance  the  effect  of  a  mag¬ 
netic  field.  A  very  recent  calculation  of 
Bdal-Monod  and  Daniel11  gives  a  complete 
analysis  of  the  scaling  factors  in  the  presence 
of  spin  fluctuations  (finite  temperature  cor¬ 
rections  to  the  T  -  0  results  of  Wohlfarth  and 
Rhodes12).  These  can  be  used  to  obtain  inter¬ 
esting  field-induced  alterations  of  the  thermo¬ 
dynamics  of  exchange  enhanced  systems. 

Using  a  method  previously  described,13  the 
field  scaling  factors  are  found  to  be: 

ST/Tf  ,  (7) 

and 

S3/2H/Tp  i  S1/2H/Ts  f  .  (8) 

Here,  S  is  the  Stoner  factor, 

S  -  (1  -  I)'1 

where 

I  -  IN(Ej.) 

We  have  taken  a  short-range  approximation  for 
the  electron-electron  lnterectlon,  I,  and  N(EF) 
is  the  density  of  states  at  the  Fermi  energy. 

In  units  of  the  square  of  the  magnetic  moment, 
the  Pauli  susceptibility  (bare)  equals 

xp.uli  ’  2N(V  • 

We  have  used  in  Eq.  (8)  the  usual  expression 
for  the  spin-fluctuation  temperature 


Beal-Monod  and  Daniel  give  the  following 
expression  for  the  magnetization  of  an 


exchange  enhanced  Fermi  gas  (parabolic  band) 
in  the  temperature  and  field  regimes  T  <<  T  ,  , 
1/2  s*‘ ' 

H  «  TS't'HS)L,i: 

M(T*H)  -  S  *PaullH  f1  -  V2(T/Tf)2 
-  80S3(H/Tf)2 

+  (8Z  +  4o1S0)  S2(T2/Tf2)  S3(H2/Tf2) 

+  •••  |  .  (9) 

Equation  (9)  is  written  in  such  a  way  that  the 
scaling  relationships  are  made  explicit.  For 
the  parabolic  band,  the  coefficients  in  Eq.  (9) 
equal 

80  -  1/6 
and  for  H  <<  T, 

o1  «  it2/6  ,  8,  »  23  it2/ (24) 2 

while  for  T  «  H, 

nz  -  tt 2 /4  ,  g  -  27  tt2/  (24)2 

Some  discussion  is  in  order.  The  T  •  0  velue 
for  S,  So,  is  the  value  computed  in  the 
Stoner-Wohlf arth  theory.12  The  temperature 
dependent  contributions  to  M(T,H)  diverge  with 
increesing  S.  Fluctuations  greatly  enhance 
finite  temperature  corrections  to  M(T,H) .  More 
detailed  discussions  can  be  found  in  the  origi¬ 
nal  and  complete  work  of  Ref.  11.  For  more 
complicated  bands  (but  still  isotropic)  O]  and 
Bq  are  known.  12’1‘‘  BBal-Monod  and  Daniel  sug¬ 
gest  that  O],  Bo,  and  Bj  will  all  have  the  same 
sign  for  arbitrary  bend  structure. 


Given  Eq.  (9),  Beal-Monod  and  Daniel11  go  on  to 
calculate  the  field  dependence  of  the  specific 
heat  coefficient  using  a  Maxwell  relation.  For 
a  parabolic  band  they  find 


y(H)  -  y (0)  -  -  xp.ull 


S3/2H 


1  - 

'el 

2^+  260 

[s3/2h 

2 

T„ 

1 

r 

. 

(10) 


Bial-Monod  points  out1'  that  the  first  term  in 
Eq.  (10)  arises  from  the  curvature  of  the  zero 
field  susceptibility  at  T  -  0.  At  high  fields, 
the  second  term  can  contribute  significantly, 
perhaps  even  reversing  the  sign  of  the  field 
dependence  of  v(H),  though  of  course  higher 
order  terms  must  also  be  Included. 


Comparison  with  experiment  is  becoming  possible 
with  the  advent  of  high  field  measurements  on 
exchanged  enhanced  materials.  The  case  of  UAlj 
will  be  analyzed  at  this  Symposium  by  F.  R. 
da  Boer  at  al.,  and  seems  to  show  the  same 
qualitative  behavior  at  predicted  by  Eq.  (9). 


6 


R.  Orbach 


B6al-Monod  and  Danial11  hava  analyzed  In  some 
decail  the  caae  of  liquid  Ha3  vhera  no  band 
structure  af facta  ara  axpectad.  At  tha  malting 
praaaura,  S  >  20,  and  tha  dapartura  from 
linaarlty  of  M(T,H)  la  predicted  to  be  of  the 
order  of  2X  at  10  Teala.  In  confined  geome- 
triea,  however,  S  can  be  meda  aa  high  aa  60, 
giving  25X  eff acta  under  the  same  conditions ■ 

The  case  of  Pd  continues  to  be  difficult  to 
sort  out.  The  susceptibility  at  H  ■  0  incraasea 
with  increasing  temperature.  This  implies, 
from  Eq.  (9),  that  aj  is  negative.  This  then 
leads  one  to  expect  an  increase  of  the  specific 
heat  coefficient  y(H)  with  increasing  H  using 
Eq.  (10).  This  seems  to  be  at  variance  with 
experiment,  though  more  recent  studies  do  seem 
to  exhibit  smaller  dacreasea  of  y(H)  with  in¬ 
creasing  H  chan  before. 14 

The  utility  of  Eqs.  (9)  and  (10)  lies  with 
their  relationship  to  one  another.  The  two  ef¬ 
fects  are  noc  independent,  and  as  seen  in  the 
example  of  Pd  metal,  there  is  a  consistency 
requirement.  One  cannot  simply  introduca  arbi¬ 
trary  coefficients  for  the  temperature  and 
fiald  dependences  of  the  susceptibility  and 
specific  heat.  Rather,  tha  various  behaviors 
ara  linked  through  a  sec  of  known  relationships. 
This  should  greatly  assist  experimental 
analysis,  and  may  serve  as  an  indication  of 
unwanted  impurities  present  when  the  consis¬ 
tency  relations  are  not  satisfied. 

4.  P-STATE  PAIRING  IN  THIN  FILM  OR  LAYERED 
SUPERCONDUCTORS 

Use  of  high  magnetic  fields  to  achieve  the  p- 
wave  condensation  state  for  superconductors  has 
been  re-examined  recently  by  Klemm  and 
Scharnberg. 1 7  There  are  many  problems  associ¬ 
ated  with  observation  of  this  state,  tha  most 
serious  perhaps  being  that  ordinary  impurity 
scattering  acts  as  a  pair-breaker  for  p-wava18 
condensation,  in  atrong  contrast  to  the  s-wava 
case  where  the  Anderson  theorem  shows  that  the 
critical  temperature  is  essentially  unaffected. 
Added  to  this  difficulty  is  the  expectation 
that  the  p-wave  transition  temperature,  T^ 

(l.e.,  triplet),  is  expected  to  be  much 
smaller  than  T.  (l.e.,  singlet)  for  s-wave  con¬ 
densation.  One  argument  favoring  triplet  pairing 
in  a  magnetic  field  is  that  singlet  pairing 
(clean,  type  II),  is  limited  by  Pauli  pair 
breaking,*^  but  triplet  pairing  la  not.  Though 
this  be  true,  Sharnberg  and  Kleom  have  recently 
shown20  that  orbital  pair  breaking  in  tha 
preaence  of  a  field  limits  triplet  condensation 
in  nearly  the  same  manner  as  for  singlet  conden¬ 
sation.  As  a  consequence,  for  bulk  materials, 
unless  Tj  Is  very  close  to  T$  (unlikely)  orbital 

pair  breaking  would  prevent  the  upper  critical 
field  H  ,  for  triplet  pairing  from  exceeding  the 
Pauli  limit. 

Tha  Issue,  then,  Is  how  to  achieve  a  condition 
where  orbital  pslr  breaking  can  be  suppressed, 


in  the  presence  of  large  enough  magnetic 
field  to  quench  singlet  superconductivity. 

Tha  idea  of  using  reduced  dimensionality  for 
this  purpose  was  first  put  forward  by  Efetov 
and  Larkin . 2 1  They  examined  layered  compounds 
[specifically,  TaS2(pyridine),  where  a  mag¬ 
netic  field  parallel  to  the  layers  is  able 
to  Induce  vortices  between  the  metallic  layers, 
and  hence  "decouple"  the  layers  from  one 
another.  This  decoupling  leads  to  an  enhance¬ 
ment  of  the  upper  critical  field  over  the  bulk 
value  by  (roughly) 

•v  itr  CT(0)/d2  (11) 

where  l[T  is  the  transport  mean  free  path,  (T(0) 

is  the  zero-temperature  triplet  coherence 
length,  and  d  the  thickness  of  the  metallic 
layer.  This  enhancement  can  be  sufficient  to 
exceed  the  Pauli  limit,  allowing  for  quenching 
of  the  singlet  state.  Unfortunately,  this  par¬ 
ticular  system  has  been  shown  to  be  very  dirty, 
leading  Klemm  and  Scharnberg  to  question  this 
explanation  for  the  observed  very  large  Hc2,ir 

In  their  paper,17  Klemm  and  Scharnberg  have 
analyzed  the  nonlocal  Gor'kov  gap  equations  for 
triplet  pairing  in  the  presence  of  magnetic 
fields  under  conditions  of  reduced  dimensional¬ 
ity. 

The  conditions  are  stringent  (specular  surface 
scattering,  clean  thin  films),  but  their  results 
can  conveniently  be  summarized  by  their  Fig.  3: 


Fig.  2.  Schematic  plot  of  H  ,  , . (T)  for  a  thin 
film.  The  shaded  region  is  * '  the  regime  of 
p-vsve  superconductivity.  Hp  is  the  Pauli 
limiting  critical  field  for  s-wave  condensa¬ 
tion.  A  similar,  though  more  complicated,  curve 
is  exhibited  by  Ref.  17  for  layered  compounds. 

Klemm  and  Scharnberg  go  on  to  suggest  physical 
syatams  which  might  exhibit  triplet  auperconduc- 
tlvlty  in  high  magnetic  fields.  They  suggest 
clesved  thin  films  of  NbSe2,  or  intercalates  of 

che  same  system. 

On*  esn  legitimately  ask,  beyond  the  structure 
of  Fig.  2,  what  algnatura  triplet  superconduc¬ 
tivity  will  give  which  is  unique  to  the  p-wave 
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paired  teat*.  Buchholtz  and  Zwlcknagl  hava 
recently  pointed  out  that  spin-polarized  tun¬ 
neling  (pioneered  by  Tedrov  and  Meiarvay)  in  a 
magnetic  field  appears  to  be  one  of  the  few 
methods  of  qualitatively  distinguishing  its 
behavior  from  that  of  tha  singlet  paired 
state.22  Another  difference  would  be  the  ef¬ 
fect  of  magnetic  proximity.  Whereas  a  ferro¬ 
magnetic  layer  in  proximity  to  a  singlet  paired 
superconductor  destroys  the  condensation  over 
a  coherence  length,  quite  the  opposite  is 
expected  to  happen  for  triplet  pairing.  In¬ 
deed,  it  might  be  feasible  to  enhance  Che 
effect  of  an  external  field  by  using  the  proxi¬ 
mity  effect  in  conjunction  with  superposed  fer- 
romagneclc  layers. 

5.  CONCLUSION 

This  paper  has  not  attempted  a  general  survey 
of  the  physics  that  can  be  done  in  high  magnetic 
fields.  Rather,  it  has  focused  on  only  three 
areas,  in  addition  to  listing  those  other  topics 
which  will  be  covered  in  depth  by  other  invited 
speakers  at  this  Sympoalw. 

The  author  suggests  that  the  main  thrust  of  the 

three  areas  he  has  discussed  are: 

1)  High  frequency  (infra-red?)  electron  spin 
resonance  relaxation  time  meaeuraments,  as  a 

function  of  frequency  in  the  direct  process 
regime  (low  temperatures,  high  fields),  to 
probe  the  "fracton"  density  of  states.  Macro- 
molecules,  cross-linked  polymare,  gels,  and 
other  "open"  self-similar  structuraa  would  be 
good  candidates  for  significant  departures  from 
tha  Debye  density  of  states  for  lattice  vibra¬ 
tions. 

2)  High  magnetic  field  studies  of  magnetization 
and  specific  heat  of  exchange  enhanced  metals, 
to  determine  the  thermodynamic  parameters  a  and 
£.  As  suggested  by  B4al-Monod  and  Daniel,  use 
of  He3  in  restricted  geometries  may  be  one  tech¬ 
nique  for  seeing  the  nonlinear  magnetic  field 
properties  of  the  Fermi  gas.  Another  may  si^>ly 
be  use  of  higher  fields  and  lower  temperatures. 

3)  Triplet  (p-wave)  superconductivity— the 
existence  of  the  state  and  its  properties. 

Recent  work  of  Klaam  and  Scharnberg  have  given 
tha  limits  on  the  magnetic  field  range  (for 
clean  syatems,  but  still  type- II)  in  which  the 
a-wave  state  is  quenched,  but  the  magnetic  field 
is  below  the  upper  critical  field  for  p-wave 
pairing.  The  uae  of  restricted  geometry  is 
essential  in  order  to  avoid  orbital  field- 
induced  pair-breaking  for  the  p-wave  state. 

Thin  films  and  layered  compounds  are  suggested 
as  excellent  candidates. 
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